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SUMMARY 


w  This  report  is  a  reference  manual  for  the  third  version  of  the  USA-STAGS  Code 
that  calculates  the  nonlinear  transient  response  of  a  totally  or  partially  submerged 
structure  to  a  spherical  shock  wave  of  arbitrary  pressure  profile  and  source  location. 
USA-STAGS  is  the  result  of  Interfacing  USA  (Underwater  Shock  Analysis)  and  STAGS  (Stress 
Analysis  of  General  Shells).  The  fluid  Is  assumed  to  be  an  Infinite  acoustic  medium 
whose  response  to  motions  of  the  structure  is  described  by  either  of  the  Doubly 
Asymptotic  Approximations,  DAA^  or  DAA^.  __ 
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SECTION  I 


INTRODUCTION 

This  report  is  a  reference  manual  for  the  third  version  of  the  USA-STAGS  Code  that 
calculates  the  nonlinear  transient  response  of  a  totally  or  partially  submerged  struc¬ 
ture  to  a  spherical  shock  wave  of  arbitrary  pressure  profile  and  source  location.  This 
constitutes  a  revision  of  [1],  the  original  manual  that  describes  the  interfacing  of  USA 
(Underwater  Shock  Analysis)  [2],  and  STAGS  (STress  Analysis  of  General  Shells)  [3]. 

The  USA  Code  computes  the  transient  response  of  a  totally  or  partially  submerged, 
shock-w.tve-^xcited  elastic  structure  using  a  linear  (or  linearised)  stiffness  matrix 
that  has  been  constructed  by  a  finite-element  or  finite-difference  structural  analyser 
of  choice.  STAGS  is  a  general-purpose  nonlinear  finite-element,  structural  code  that  is 
very  efficient  for  the  analysis  of  inelastic  collapse  of  stiffened  shell  structures. 
USA-STAGS  employs  the  staggered  solution  strategy  upon  which  USA  itself  is  con¬ 
structed.  In  this  procedure  the  structural  response  equations  and  the  fluid  response 
equations  are  solved  separately,  step-by-step,  through  extrapolation  of  terms  which 
couple  the  two  systems  thus  preserving  the  modularity  of  each  code.  This  technique  has 
been  shown  to  be  unconditionally  stable  with  respect  to  time  step  siae  in  the  linear 
response  problem  [4]. 

The  fluid  equations  treated  by  USA-STAGS  apply  to  an  infinite  or  semi-infinite 
acoustic  medium  and  either  the  Doubly  Asymptotic  Approximation  (DAAj,)  or  the  Improved 
Doubly  Asymptotic  Approximation  (DAA^ )  [5,  6]  is  used  to  determine  the  fluid  pressure  on 
the  "wet"  surface  of  the  structure. 

The  following  additions,  changes  snd  enhancements  have  been  incorporated  in  the 
USA-STAGS  Code  since  [1]  was  issued: 

1)  STAGS  C-l  has  become  the  standard  version  for  underwater  shock  analysis. 

2)  An  imaging  formulation  that  allows  the  analysis  of  problems  involving  partially 
submerged  structures  and  structures  totally  submerged  near  the  free  surface  ’.ias 
been  implemented. 

3)  A  technique  has  been  developed  to  include  the  effects  of  bulk  cavitation  in  the 
incident  wave  excitation  that  provides  a  complete  description  of  fluid  particle 
velocity  consistent  with  the  occurrence  of  surface  cutoff. 

4)  Surface-of-revolution  fluid  elements  can  now  treat  both  beam-like  and  bar-like 
response. 

5)  The  Improved  Doubly  Asymptotic  Approximation  (DAA£)  is  now  operational. 


6)  The  fluid  equation  system  is  processed  in  an  out-of-core  mode  to  allow 
substantial  increases  in  the  allowable  number  of  fluid  elements. 

7)  A  capability  to  do  linear  analysis  has  been  built  into  the  USA-STAGS 
Code. 

8)  An  improvement  has  been  added  to  the  integration  time  step  control  to  allow 
STAGS  to  change  the  step  siae  if  required  for  convergence. 

9)  A  hydrostatic  prestress  analysis  of  the  structure  can  be  made  tc  provide  non- 
zero  initial  values  for  nonlinear  underwater  shock  analysis. 

10)  A  capability  to  conduct  a  collapse  analysis  of  the  structure  upon  completion 
of  the  transient  analysis  has  been  added. 

11)  USA-STAGS  is  now  operational  on  the  VAX  11/780  virtual  memory  system  machine. 

1.1  STAGGERED  SOLUTION  PROCEDURE 

USA-STAGS  is  constructed  around  the  staggered  solution  procedure  which  avoids  the 
simultaneous  solution  of  the  structural  response  equations  and  the  fluid  response  equa¬ 
tions.  Such  a  simultaneous  solution  would  be  costly  from  two  viewpoints.  First,  a 
special-purpose  version  of  STAGS  would  have  b  m  required.  Second,  the  computation 
time  would  have  been  substantially  increased  over  that  required  for  the  transient 
analysis  of  a  dry  structure  due  to  the  widespread  coupling  induced  by  the  fluid. 

On  the  other  hand,  the  staggered  solution  procedure  involves  a  response  extrapola¬ 
tion  at  each  time  step,  which  usually  leads  to  numerical  instability  for  time  incre¬ 
ments  exceeding  a  critical  value.  Because  this  critical  value  may  be  unacceptably 
small  for  some  computations,  the  fluid  equations  used  in  USA  have  been  modified  in  such 
a  way  that  unconditional  stability  has  been  achieved  for  the  linear  elastic  structural 
response  problem  [4], 

At  this  time  an  investigation  of  time  step  stability  of  the  coupled  USA-STAGS  sys¬ 
tem  has  not  been  performed.  Stability  will  depend  not  only  upon  the  specific  types  of 
nonlinearity  inherent  in  the  STAGS  equations  but  also  upon  which  of  the  family  of  STAGS 
integration  algorithms  is  used  in  any  particular  case.  It  is  likely,  however,  that 
unconditional  stability  will  exist  for  some  range  of  mildly  nonlinear  behavior. 

1.2  DOUBLY  ASYMPTOTIC  APPROXIMATIONS 

The  DAA^is  embodied  in  a  linear  first  order  matrix  ordinary  differential  equation 
that  defines  the  fluid  pressure  due  to  the  scattered  wave  at  the  fluid-structure  bound¬ 
ary.  The  principal  advantage  of  the  DAA^then  is  that  it  models  the  infinite  acoustic 
medium  surrounding  the  structure  as  a  membrane  covering  the  wet  surface  of  the  struc- 
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tare.  Hence  fluid  motion  Is  described  merely  In  terms  of  wet-surface  response  varia¬ 
bles,  which  are  then  linked  by  compatibility  relations  to  the  structural  response  varia¬ 
bles. 

The  principal  disadvantage  of  the  DAAj  is  that  it  constitutes  an  approximation  to 
the  "exact"  boundary-element  representation  of  the  surrounding  medium  [7,8].  The  DAAj 
does  approach  exactness  for  both  high-frequency  (early-time)  and  low-frequency  (late¬ 
time)  structural  motions,  however,  and  effects  a  smooth  transition  between  the  two 
asymptotes.  In  addition,  it  has  exhibited  satisfactory  accuracy  in  a  variety  of  check 
calculations  [6,7,9]. 

The  DAA2  is  an  improved  approximation  that  Is  based  upon  the  DAAj ;  however,  it  can 
describe  the  intermediate  frequency  range  better  than  the  DAAj  because  it  is  a  second 
order  matrix  differential  equation.  Computationally  it  has  been  used  to  study  the  res¬ 
ponse  of  the  infinite  cylindrical  shell  under  a  plane  wave  step  lording  in  which  signi¬ 
ficant  improvements  in  accuracy  have  been  noted  [2],  Analytical  studies  of  spherical 
shell  response  [6]  also  show  such  marked  improvement. 

Although  the  DAAj  and  DAA2  have  been  shown  to  be  suitable  for  engineering  analysis 
of  underwater  shock  problems,  the  physics  of  such  problems  still  include  features  of 
great  complexity  not  currently  treated.  For  example,  local  fluid  cavitation  at  the 
fluid-structure  interface  may  substantially  alter  structural  response  for  incident  shock 
waves  of  short  duration.  Effective  means  of  including  hull  cavitation  phenomena  in  the 
numerical  solution  of  underwater  shock  problems  are  still  being  sought. 


1-3 


U3  USA  FEATURES 


A  number  of  special  features  are  incorporated  In  the  USA  code.  First,  a  capability 
has  been  provided  to  handle  a  fluid  mesh  on  the  wet  surface  that  is  not  coincident  with 
the  surface  mesh  for  the  structural  models  This  permits,  for  example,  the  use  of  a 
refined  structural  mesh  in  a  region  of  high  stress  gradients,  even  though  a  relatively 
coarse  mesh  is  retained  for  the  fluid. 

Second,  options  for  variable-increment  time  integration  and  computation  restart  are 
furnished.  The  former  allows  the  use  of  small  time  increments  during  periods  where  the 
response  is  expected  to  be  varying  rapidly  in  time,  and  the  use  of  large  time  increments 
tor  periods  characterised  by  a  slowly  varying  response.  The  latter  permits  the  division 
of  a  response  computation  into  segments,  so  that  the  analyst  may  examine  the  results  at 
selected  points  along  the  way.  Such  examination  is  facilitated  by  the  use  of  the 
"printer-plot"  routine  that  augments  the  usual  printout  data  with  response  plots  "drawn" 
by  the  printer. 

Third,  the  code  incorporates  fluid  boundary  elements  for  both  general  and  body-of- 
revolution  wet-surface  geometries  [10],  This  feature  is  especially  useful  for  compart- 
ment-by-compartment  analysis  of  a  submarine.  Such  an  analysis  utilizes  a  general- 
structure  discrete-element  model  of  a  particular  compartment  of  interest,  with  the  re¬ 
mainder  of  the  submarine  modeled  as  a  bar/beam.  Hence  a  detailed  analysis  of  an  entire 
submarine  may  be  performed  with  several  discrete-element  models  of  moderate  size, 
avoiding  the  use  of  a  single  gigantic  model. 

Fourth,  the  analyst  can  use  either  the  DAAj  or  DAA2  through  input  of  a  single 
scalar  parameter  to  take  advantage  of  the  enhanced  accuracy  demonstrated  by  the  latter 
in  computations  for  idealized  geometries. 

Fifth,  out-of-core  processing  for  both  the  structural  and  fluid  equation  systems 
frees  the  user  from  concern  over  core  limitations  on  the  number  of  structural  and  fluid 
elements  in  his  model. 

Sixth,  free-field  shock  wave  input  to  the  structure  is  associated  with  a  spherical 
wave  for  both  submerged  and  partially  submerged  structures  and  can  be  input  for  use  with 
linear  interpolation  or  cubic  spline  fitting  routines.  Pressure  hisories  for  exponen¬ 
tially  decaying  incident  waves  can  be  automatically  generated.  Fluid  pressure  and  par¬ 
ticle  velocity  histories  corresponding  to  the  input  shock  are  displayed  for  the  user 
with  the  "printer-plot"  software  mentioned  above. 

Seventh,  the  effects  of  bulk  cavitation  on  the  free  field  pressure  history  are 
approximately  treated. 
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1.4  STAGS  FEATURES 


STAGS  is  primarily  intended  for  the  nonlinear  analysis  of  shell  structures; 
therefore,  great  emphasis  has  been  placed  on  computational  efficiency.  The  capabili¬ 
ties  of  STAGS  include  the  following  processors:  1)  linear  stress  analysis;  2)  geomet¬ 
rically  nonlinear  elastic  stress  analysis;  3)  inelastic  stress  analysis,  geometrically 
linear  or  nonlinear;  4)  bifurcation  buckling  analysis  with  linear  or  geometrically  non¬ 
linear  prestress  (elastic);  5)  small  vibration  analysis  with  prestress  based  on  linear 
or  geometrically  nonlinear  analysis  (elastic);  and  6)  transient  response  analysis, 
linear  or  geometrically  nonlinear,  elastic  or  inelastic.  In  a  STAGS  analysis,  the 
structure  can  be  defined  as  composed  of  a  number  of  shell  branches.  Linear  springs 
(axial  or  torsional)  and  nonlinear  general  beam  elements  also  can  be  included.  If  the 
beams  are  defined  as  stiffeners  attached  to  the  shell  along  a  line  on  its  surface, 
appropriate  displacement  constraints  will  be  enforced  automatically.  Loading  is  either 
mechanical  (forces  or  forced  displacements)  or  thermal.  Any  boundary  conditions  or 
other  displacement  constraints  are  permitted  as  long  as  they  are  linear.  Displacement- 
dependent  hydrostatic  loads  may  be  used.  For  transient  problems,  velocity-dependent 
damping  may  be  included  in  specified  regions.  General  damping  proportional  to  mass  or 
linear  stiffness  may  also  be  selected. 

A  number  of  quadrilateral  and  triangular  finite-elements  are  available  as  options 
for  discretization  of  the  shell  branches.  The  discretization  procedures  available  in¬ 
clude  a  flat  element  in  which  lateral  displacements  are  bicubic  and  in-plane  displace¬ 
ments  at  least  quadratic.  At  each  corner  node,  there  are  six  or  seven  degrees  of  free¬ 
dom:  three  displacement  components,  three  rotation  components,  and  optionally  the  shear 
strain.  In  addition,  displacements  along  the  element  boundaries  may  be  included  as 
degrees  of  freedom  at  midpoint  nodes,  thus  the  maximum  number  of  degrees  of  freedom 
per  element  is  32.  Refined  quadrilateral  elements  constructed  from  sub-triangles  may 
also  be  used.  These  quadrilaterals  need  not  be  flat. 

Bifurcation  buckling  and  small  vibration  analyses  lead  to  linear  eigenvalue  prob¬ 
lems.  In  STAGS,  these  are  solved  through  the  generation  of  invariant  subspaces  by 
simultaneous  power  iteration  [11].  Nonlinear  static  analysis  results  in  nonlinear 
algebraic  equations.  For  solution  of  such  equations,  the  STAGS  user  has  the  options 
of  using  either  the  modified  or  the  regular  Newton  method  [12].  For  transient  analysis 
numerical  integration  the  following  options  are  available:  1)  the  central  difference 
scheme;  2)  trapezoidal  rule;  3)  Gear's  two-  and  three-order  methods;  or  4)  Park’s 
method.  A  discussion  of  these  methods  is  available  in  Ref.  [13]. 
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1.5  SPECIAL  NON-FEATURES 


Some  features  of  modest  complexity  have  yet  to  be  Incorporated  into  the  USA-STAGS 
Code.  First,  an  option  for  automatic  time-step  integration  would  free  the  analyst  from 
having  to  select  integration  time  increments  in  accordance  with  his  expectations 
regarding  response  behavior.  Second,  the  ability  to  handle  very  large  problems  would  be 
useful  in  those  cases  where  structural  segmentation  is  not  possible.  Third,  a  capa¬ 
bility  to  treat  banded  structural  mass  and  damping  matrices  would  be  helpful,  in  order 
to  accommodate  structural  elements  that  produce  such  matrices.  Fourth,  a  means  to 
handle  the  matrices  produced  by  acoustic  elements  based  on  a  pressure  formulation 
[14,15]  is  desirable,  in  that  such  elements  permit  highly  efficient  modeling  of  internal 
fluid  volumes. 

An  important  feature  of  greater  complexity  has  yet  to  be  introduced  into  the 
code.  This  is  a  treatment  of  hull  cavitation,  which  may  substantially  affect  structural 
response  for  incident  shock  waves  of  short  duration.  The  introduction  of  this  feature 
requires  the  accurate  treatment  of  highly  nonlinear  phenomena,  and  presents  a  challen¬ 
ging  task  for  future  work. 
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SECTION  II 


FLUID-STRUCTURE  INTERACTION  PROBLEM  FORMULATION 
AND  SOLUTION  PROCEDURE 

This  section  describes  the  theoretical  foundation  of  the  USA-STAGS  Code;  and  la  con¬ 
structed  as  an  overview,  with  coverage  of  details  left  to  referenced  papers  and  reports. 

2.1  STRUCTURAL  RESPONSE  EQUATION 

The  discretized  differential  equation  system  for  the  dynamic  response  of  a  nonlinear 
structure  can  be  expressed  in  the  form 


Mx  +  Cx  +  Kx-f 
—s  —  ~s  —  ~s  —  — 


(2.1) 


where  x  is  the  structural  displacement  vector,  M  ,  C  and  K  are  the  nonlinear  structural 
—  ~ S  ~S  '—8 

mass,  damping  and  stiffness  operators,  respectively,  f  is  the  external  force  vector,  and 
a  dot  denotes  a  temporal  derivative.  Generally,  M  ,  C  and  K  are  highly  banded,  symmetric 
operators  of  large  order;  at  present  the  USA-STAGS  Code  considers  to  be  diagonal  and  (Cg 
to  be  a  linear  combination  of  Ms  and  the  linear  portion  of  Kg. 


For  excitation  of  a  submerged  structure  by  an  acoustic  wave,  f  is  given  by 

f  *  -£  Af  (El  +  £g)  +  iD  (2'2) 

where  gj  and  gg  are  nodal  pressure  vectors  for  the  wet-surface  fluid  mesh  pertaining  to  the 
(known)  incident  wave  and  the  (unknown)  scattered  wave,  respectively,  fp  is  the  dry-struc¬ 
ture  applied  force  vector,  A^  is  the  diagonal  area  matrix  associated  with  elements  in  the 
fluid  mesh,  and  £  is  the  transformation  matrix  that  relates  the  structural  and  fluid  nodal 
surface  forces.  More  will  be  said  about  G,  in  the  next  subsection. 

2.2  DAAX  EQUATION 

The  Doubly  Asymptotic  Approximation  may  be  written  [5,6] 


Mf  Pg  +  pc  Af  Pg  -  pc  Mf  us 


(2.3) 


where  ug  is  the  vector  of  scattered-wave  fluid-particle  velocities  normal  to  the 
structure's  wet  surface,  p  and  c  are  the  density  and  sound  velocity  of  the  fluid, 
respectively,  and  is  the  symmetric  fluid  mass  matrix  for  the  wet-surface  fluid  mesh. 
This  matrix  is  produced  by  a  boundary-element  treatment  of  Laplace's  equation  for  the 
irrotational  flow  generated  in  an  infinite,  inviscid,  incompr es sibl e  fluid 
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by  motion*  of  th*  structure's  wet  surface;  It  is  fully  populated  with  non-zero  matrix  ele¬ 
ments.  When  transformed  into  structural  coordinates,  the  fluid  mass  matrix  yields  the  added 
mass  matrix,  which,  when  combined  with  the  structural  mass  matrix,  yields  the  virtual  mass 
matrix  for  motions  of  a  structure  submerged  in  an  incompressible  fluid  (lb). 


Ac  mentioned  in  Section  I,  the  approximate  relation  (2. 3)  is  called  "doubly  asymptotic" 
because  it  approaches  exactness  in  both  the  high-frequency  (early-time)  and  low-frequency 
(late-time)  limits.  For  high-frequency  motions,  |£g|  >>  jp^l  ,  so  that  (2.3)  approaches  the 
relation  £g  •  pcu^  ,  which  is  the  correct  limit  for  short  acoustic  wavelengths.  For  low- 
frequency  motions,  l^l  «  |£g|  ,  so  that  (2. 3)  approaches  the  incompressible-flow  relation 
A^  Ug  ,  which  is  the  correct  limit  for  long  acoustic  wavelengths. 

For  excitation  by  an  incident  acoustic  wave,  Ug  is  related  to  structural  response  by 
the  kinematic  compatibility  relation 


T. 

64 


Hi  +% 


(2.4) 


where  the  superscript  "T"  denotes  matrix  transposition.  Equation  p.Q  expresses  the  con¬ 
straint  that  normal  fluid-particle  velocity  match  normal  structural  velocity  on  the  wet 
surface  of  the  structure.  The  fact  that  the  transformation  matrix  relating  those  velo- 
cities  is  G  follows  from  the  invariance  of  virtual  work  with  respect  to  either  of  the 
wet  surface  coordinate  systems.  Generally,  G  is  a  rectangular  matrix  whose  height 
greatly  exceeds  its  width,  inasmuch  as  the  number  of  structural  DOF  usually  exceeds  con¬ 
siderably  the  number  of  fluid  DOF. 


2 . 3  INTERACTION  EQUATIONS 


The  introduction  of  (2.2)  into  (2.1)  and  (2.4)  into  (2.3)  yields  the  interaction  equa¬ 
tions 


Ms  I  +  Cs  x  +  Ks  x  -  -  G  Af  <El  +  Es) 

Mf  Es  +  p  C  ~f  Ps  "  P  c  ~f  (~TH  "  V 


(2.5)* 


These  equations  may  be  solved  simultaneously  at  each  time  step  by  the  transfer  of  -  G  A^,  p^ 
T 

and  p  c  G  x  to  the  left  sides  of  their  respective  equations.  Such  a  procedure  is 
~f  ~  — 

exceedingly  expensive,  however,  because  of  the  large  connectivity  of  the  coefficient  matrix 
involved.  As  mentioned  in  Section  I,  efficient  computation  is  possible  through  the  appli¬ 
cation  of  a  staggered  solution  procedure  that  is  unconditionally  stable  with  respect  to 
the  choice  of  time  increment,  at  least  in  the  linear  case. 

The  simplest  implementation  of  the  staggered  solution  procedure  recommended  in  [4]  may 
be  effected  as  follows.  M„  is  taken  to  be  diagonal  and,  to  allow  for  the  possibility  that 
M  may  have  zero  entries  for  rotational  DOF,  £  is  constructed  such  that  only  the  transla- 
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Clonal  DOF  for  the  structure  couple  with  the  fluid  DOF  [see  (2.4)];  then  the  first  of  (2.5) 
nay  be  partitioned  to  obtain  x,  which  Bay  then  be  Introduced  Into  the  second  of  (2.5). 
Premultiplication  of  the  resulting  equation  by  oc  Af  ~f  then  yields 


n  +  <Sn+  £s)  *s 


-Af  GT  x  +  Kfl  x) 


“(Ba  i-i  +  Af  %) 


(2.6) 


where  D^**  A^  and  Dg  ~  GA  M~A  G  A^  are  symmetric,  and  where  Ms"  is  a  diagonal 

matrix  with  each  nonzero  element  given  as  the  reciprocal  of  the  corresponding  nonzero  element 
of  Mg  and  each  zero  element  mirroring  the  corresponding  zero  element  of  Ms  .  The  first 
of  (2.5)  and  (2.6)  are  herein  termed  "the  augmented  interaction  equations". 

2.4  SPHERICAL  INCIDENT  WAVE 

Each  element  of  the  vectors  and  Uj  for  a  spherical  incident  wave  are  given  by 


T  -1 

1*  Vi  * 


S  /  Ri ~S  \ 

Pn(t)  "  T{  p!  -  C  ) 

uli(t)  ■  [  pc  +  pRj  pIi^^  Yi 


(2.7) 


whe-re  S  is  the  "charge  standoff",  i.e.,  the  distance  between  the  origin  of  the  incident 
spherical  wave  and  the  nearest  point  on  the  structure's  wet  surface,  R^  is  the  distance 
from  the  origin  of  the  incident  spherical  wave  to  the  ith  fluid  node  on  the  wet  surface, 
y^  is  the  cosine  of  the  angle  between  the  vector  corresponding  to  R^  and  the  wet-surface 
normal  at  the  i£li  fluid  node,  and  Pj(t)  is  the  incident-wave  pressure-profile  defined  at 
*  S  .  For  3  shock  wave,  p^(t)  *s  discontinuous  at  t  ■  0  and  the  u^(t)  contain 
singularities. 

In  order  to  remove  shock-wave  singularities  from  u^  in  (2.6),  a  modified  pressure 


vector  is  defined  as 


Em  •  1&I  +  *s 


(2.8) 


where  r  is  a  diagonal  matrix  with  direction-cosine  elements  .  The  introduction  of 

(2.8)  into  (2.6)  and  the  first  of  (2.5),  followed  by  utilization  of  the  second  of  (2.7)  then 
yields  the  modified,  augmented,  interaction  equations 


Us  2  +  £g  5  +  h  5  "  -  £  Af  [Em  +  (I  -  r)  Pj] 


(2.9) 


7E  Af  Em  +  (£fi+  Ss  )  em  =  -  *t  £T  ifs1  x  +  k8  x)  -  hEi 
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in  which  jl  Is  the  identity  matrix,  end 

a  ■  e.  -  <£. +  fin-  i  ^  £_1>  z  «-10> 

where  £  is  the  diagonal  matrix  formed  by  the  distances  R^.  Equations  (2.9) 
are  the  equations  solved  by  the  USA-STAGS  to  determine  the  structural  responses 
x  and  x,  and  the  wet-surface  pressures  £  ■  (I_  ~  £) 

2.5  FREE  SURFACE  EFFECTS 

When  a  structure  is  partially  submerged,  or  when  a  totally  submerged  structure  lies 
near  the  free  surface  of  a  semi-infinite  fluid,  imaging  techniques  may  be  utilised  to  en¬ 
sure  that  the  total  pressure  vanishes  at  the  free  surface.  (This  implies  that  the  effects 
of  gravity  are  negligible  in  this  class  of  problems,  which  they  generally  are.)  In  this 
case,  the  interactive  system  consists  of  an  infinite  fluid  domain,  the  structure  S+,  and 
its  image  S_  (see  Figure  2-1).  The  Incident  wave  now  consists  of  a  (positive)  primary 
wave  plus  a  (negative)  image  wave,  the  latter  emanating  from  the  image  of  the  primary 
wave's  origin.  Zero  pressure  at  the  free  surface  is  therefore  maintained  if  the  motions 
of  S_  are  constrained  to  be  opposite  to  those  of  S+. 

The  kinetic  energy  T  ,  the  Rayleigh  dissipation  function  D  ,  the  potential  energy  V  , 

s  s  s 

and  the  work  potential  n  ,  for  the  structural  system  S  +  S_  are  given  by 

T.  -  *  <£  4  i.  +  £  H, 

D  ■  J  (x+  C  Xx  +  x*  C  x_) 

(2.11) 

v»  -  *  <a£  £«  4  *-> 

ns  "  -2+  |+  “  5^  t~ 

The  appropriate  constraints  are  x_ * -x+  and  f_  *  -f+,  so  that  (2.11)  become 

T  -  xl  M  x+ 

s  —s  - 

D  ■  xj  C  ii 

(2.12) 

V  «  xJk  x+ 

ns  -  -2  5+  *+ 
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Figure  2-1  Image  Technique  for  Figure  2-2  Free  Field  Ray  Diagram 

Free  Surface 


The  DAA  kinetic  energy  Tf  end  work  potential  nf  for  the  fluid  system  may  be  written 


Tf  '  i  ug  Mf  ug 

*T  i  *t 

nf  "  -Hs  h[  Eg  "  ^  2S  Mf  ts 


(2.13) 


where 


rM  M’l  Tk  0" 

Mf  "  .  Af  - 

^  M'  M  ^  0  A 


(2.14) 


and  an  asterisk  denotes  temporal  Integration.  The  submatrix  J£  accounts  for  added  mass 
coupling  between  wet-surface  elements  on  S+  and,  similarly,  between  elements  on  S_;  M' 
accounts  for  added  mass  coupling  between  elements  on  S+  and  elements  on  S  .  The  con¬ 
straints  for  the  fluid  system  are 


*  2<; 


(2.15) 


so  that  (2.13)  becomes 


T e  -  uT.  (M  -  M')  u 
f  -S+  ~  ~  -s+ 

( 

"f  ■  -2  V  I*  Es*  +  £  (S  -  S’>  is*  1 

The  application  of  Lagrange's  equation  [17]  to  (2.12)  and  (2.16)  now  yields 

M  x.  +C  i+K  x  *  f 

— s  -+  ~s  -+  s  -+  -+ 


(2.16) 


~  <5-r>  is+  +  &ES+  •  <«-h'>ss+ 


(2.17) 


Also,  (2.2)  and  (2.4)  must  be  modified  to  include  the  effects  of  both  the  incident  prim¬ 
ary  and  Image  waves.  This  gives 


-+  "  ~  *Ei+  +  Ei+  +  2s+J 


£t  *+  -  £+  +  *1+  +  HS+ 


(2.18) 


viere,  e.g.,  £I+  denotes  Incident-wave  pressure  on  S+  associated  with  the  image  wave. 
The  introduction  of  (2.18)  Into  (2.17)  then  yields  the  doubly  asymptotic  interaction 
equations  for  problems  involving  a  free  surface 


(2.19) 


~s  ?♦  +  £a  »+  *  i,  «+  •  £  *  <et+  +  i~i+  *  BSt> 

<H  -  “>  Es+  +  cc  A  fs+  "  pc  <£T  ;+  -  Hjt  -  sjt) 

A  comparison  of  (2.19)  with  (2.5)  reveals  that  the  affects  of  the  free  surface  are  embod¬ 
ied  in  the  image-wave  pressure  and  f luid-particlc-velocity  vectors,  and  in  the  modified 
added -mass  matrix. 

Finally,  augmentation  of  (2.19)  to  secure  unconditional  stability,  followed  by 
introduction  of  the  modified  pressure  [cf.  (2.8)] 

bm  "  4  4+  +4  4+  +  ps+  (2.20) 

to  remove  shock-wave  singularities,  proceeds  as  described  in  Subsections  2.3  and  2.4. 

The  modified,  augmented  interaction  equations  corresponding  to  (2.9)  for  the  infinite 
fluid  medium  are  then  readily  obtained. 

It  is  important  to  mention  at  this  point  that  the  DAA  formulation  just  described 
does  not  account  for  high-frequency  scattered  waves  from  S  that  impinge  upon  S..  For 

"*  "T" 

most  floating  structures,  such  waves  are  not  generated,  as  the  wet  surfaces  of  S  and  S 

T  "• 

usually  intersect  to  form  a  convex  surface;  they  are  generated,  however,  for  a  totally 
submerged  structure  lying  near  the  free  surface.  Even  so,  it  has  been  shown  that,  as 
far  as  structural  response  is  concerned,  the  effects  of  the  scattered  wave  are  generally 
negligible  [18].  In  other  words,  the  response  is  basically  driven  by  the  incident 
primary  and  image  waves. 

2.6  BULK  CAVITATION 

In  the  absence  of  bulk  cavitation,  the  imaging  method  serves  as  a  useful  device  to 
model  the  reflection  of  free-field  waves  from  the  fluid's  free  surface.  The  occurrence 
of  bulk  cavitation  near  the  surface,  however,  changes  that  simple  acoustic  reflection 
problem  into  a  complex  reflection-refraction  problem,  as  indicated  in  Figure  2-2.  If 
refraction  distortions  produced  by  a  relatively  thin  cavitated  region  are  not  too  severe, 
however,  bulk  cavitation  effects  will  still  appear  to  the  structure  as  emanating  from  an 
image  source. 

Experimental  records  of  free-field  pressure  histories  for  compact  charges  exhibit 
the  behavior  shown  in  Figure  2-3  [19].  The  dashed  line  denotes  the  history  produced  by 
a  negative- image  model,  while  the  horizontal  line  indicates  that  the  effect  of  bulk  cavi¬ 
tation  is  to  "cut  off"  the  pressure  at  a  cavitation  threshold.  The  approximate  treatment 
introduced  here  involves  pre-examination  of  the  image-based  free-field  pressure  at  the 


k 
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standoff  point,  i.e.,  the  point  on  the  submerged  structure  closest  to  the  charge.  When* 
ever  that  pressure  becomes  negative  to  the  extent  that  its  magnitude  exceeds  the  absolute 
ambient  pressure  at  that  depth,  a  positive  contribution  is  incorporated  into  the  negative- 
image  source  so  that  the  free-field  pressure  at  the  standoff  point  never  dips  below  the 
cavitation  threshold.  When  the  structure's  overall  dimensions  are  small  relative  to  the 
durance  from  the  structure  to  the  cavitated  region,  the  effects  of  the  positive  contri¬ 
bution  will  not  vary  appreciably  in  the  vicinity  of  the  structure. 

The  preceding  discussion  leads  to  the  following  development.  The  two— source  model  of 
Figure  2-4  yields  as  the  free-field  pressure  at  any  point  P 


Pp(t) 


irp-  (t 


(2.21) 


where  p+(t)  -  p_(t)  -  0  for  t  <  0.  At  the  standoff  point,  (2.21)  becomes 

s  -  s 

Pq^)  "  P+(t)  +  jf-  P.  (t  -  - -  )  (2.22) 


Now  p_(t)  -  -p+(t)  as  long  as  the  resulting  Pq(t)  exceeds  the  cavitation  threshold  so  that 
"cutoff"  does  not  occur;  otherwise,  pQ(t)  remains  at  the  threshold  value  -(PA  +  Y  2) 
where  pA  is  atmospheric  pressure  and  y  is  the  fluid's  weight  density.  Hence,  during  the 
"cutoff  period", 

S_  S_  -  S 

P_(t)  -  "  (P+  (t  +  - - )  +  PA  +  Y  2]  (2.23) 

The  model  just  described  fits  a  prescribed  free-field  pressure  history  at  the  stand¬ 
off  point  in  such  a  way  that  surface  cutoff  effects  appear  to  the  structure  as  emanating 
from  an  image  source.  Because  the  model  is  complete,  it  also  provides  the  free-field 
fluid-particle-velocity  information  required  for  DAA  calculations.  The  usefulness 
of  the  model  has  been  demonstrated  from  the  results  of  free-field  tests  specifically 
designed  to  produce  both  pressure  and  fluid-particle-velocity  data  l20j. 

2.7  USA-OM2  IMPLEMENTATION 

The  Improved  Doubly  Asymptotic  Approximation  jDAAjJ  can  be  written  as  J6  ] 


Mf£s+pc£f£g  +  Pc  8f  Af£g  - 
Pc  Itff(GTx-U  j)  gf  Mf(GTx-  u  j)] 
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where 


Qf  -  n  pc  *£  fc"1  (2.25) 

Ml  vector  and  matrix  quant itia*  in  tha  above  era  related  to  the  same  finite  element 
vat-aurfaca  fluid  mesh  aa  that  used  for  tha  lowest  order  DAA  (DAA^)  described  in 
Section  2.2. 

Note  that  DAA^  is  a  second-order  equation,  whereas  DAA^  is  a  first-order  equation.  In 
addition,  DAA2  includes  a  new  scalar  parameter  r)  that  appears  in  (2.25).  It  can  be 
established  from  physical  considerations  [A]  that  n  must  be  bounded  as 

0  <  n  <  1  (2.26) 

A  precise  choice  of  n  is  apparently  not  predicated  by  any  fundamental  principle. 

Hence  it  must  be  regarded  at  this  time  as  a  factor  which  may  be  adjusted  to  achieve 
optimum  accuracy  for  a  particular  problem.  In  [6],  it  is  observed  that  n  *  1  leads 
to  the  best  accuracy  for  a  spherical  shell . 


In  order  to  implement  DAA-  (2.24)  is  first  integrated  once  in  time  and  multiplied 
through  by  A£  M£  .  Equation  (2.25)  is  then  substituted  into  the  result  and  a  new 
variable,  the  scattered  pressure  integral  ,  is  defined  by 


(2.27) 


where  an  asterisk  denotes  temporal  integration.  The  result  is 

*f  +  pc  Bfii-s +  np  2  c2  5 f 2  " 

p  c  6  f  <ST  *  -  A  x)  +  n  p2  c2  d  f x  (gt  x  -  u  x) 


(2.28) 


where 


(2.29) 
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It  it  noted  that  (2.28)  it  lyutetric  and  that  D^haa  already  been  defined  following 

(2.6). 


To  avoid  ahock-wave  aingularitiea  in  d  ,  the  relation  for  a  tpherical  thock  it 

used  at 


-1 


pc  dj  -  T  (£j  +  cR  £j) 


(2.30) 


while  the  modified  pressure- integral  vector  is  defined  as 


a„-as  +  £*] 


(2.31) 


Substitution  of  (2.30,  (2.31)  into  (2.28)  then  gives 


2  2 


AfIM  +  PcDfl£M+nP  c  Df2*M  - 
pcA{  GT  jc  +  n  p2  C2  D  flGT  x  +  c  [(1  -  n)  p  5  fi  "  A  f  R  1] 


(2.32) 


+  np2c2  <Df2-}5fi  5'lj  r*i 


where  the  identity 


R'1  r  -  r  l1  (2.33) 

has  been  used  in  (2.32),  as  both  matrices  are  diagonal.  Associated  with  (2.32)  is 
the  structural  equation  of  motion 

M  x+C  x  +  K  *  -  -  G  A.  IqM  +  (I -T  )  £tJ  (2.34) 

Equations  (2.32)  and  (2.34)  define  the  DAA^-modified'interaction  equations  that  are 
solved  according  to  the  staggered  solution  strategy;  hence  an  examination  of  stability 
must  be  conducted.  It  has  been  shown  that  the  step-by-step  integration  of  (2.32)  and 
(2.34)  is  conditionally  stable;  however,  no  systematic  study  of  stability  has  yet  been 
undertaken.  In  view  of  the  fact  that  unconditional  stability  was  achieved  for  USA- 
DAA1  by  augmentation,  and  that  (2.24)  is  essentially  the  DAA^  with  a  correction  term, 
augmentation  of  (2.32)  was  carried  out  in  the  same  manner  as  that  used  for  DAA^. 
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Accordingly,  (2.34)  1*  first  solved  for  it  end  substituted  into  (2.32)  ti  give 

if  in+  ” c  « ti+  5  .>  iM  +  n '2  a,  ■ 

-  P  c  A  f  GT  M~1(C  ^  x+K  #x)  +np  2c2DflGT  * 

(2. 

-  P‘{e8-t$,+  U-n)Dfl  --Af  R'1)  r)  2X 
+  r,p2c2  <5  f  2  -  I  B  £1  S_1>  tlx 

where  D  hes  elreedy  been  defined  following  (2.6). 

* *  ■ 

Equations  (2.34)  and  (2.35)  are  the  DAAj-aodified,  augmented  interaction  equations 
that  have  been  implemented  in  the  USA  Code. 
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ORGANIZATION  OF  USA-STAGS  CODE 

Organization  of  the  USA-STAGS  system  is  shown  in  Figure  3-1.  STAGS  operates  on  two 
levels  within  this  structure;  first  as  a  preprocessor,  and  second  as  a  driven  utility 
during  the  time  integration  phase  of  the  analysis.  In  the  preprocessing  mode  STAGS 
generates  two  types  of  information.  The  first  is  data  that  are  required  later  by  the 
STAGS  time  integrator,  and  the  second  is  data  that  are  needed  for  fluid  preprocessing. 
Both  sets  of  information  are  placed  in  mass  storage  where  the  latter  data  are  retrieved 
by  the  fluid  mass  matrix  processor  FLUMAS  and  then  by  the  processor  AUGMAT.  AUGMAT 
produces  the  final  form  of  the  fluid  matrix  equations  to  be  integrated.  Data  transfer 
within  USA  is  carried  out  with  the  data  manager  DMGASP,  while  out-of-core  processing  in 
FLUMAS  is  accomplished  by  the  virtual  memory  simulator  VMSYST.  Out-of-core  processing 
of  the  fluid  equation  system  during  time  integration  is  carried  out  by  the  SKYPUL  pro¬ 
cessor  which  is  not  shown  in  3-1. 

The  time  integration  processor  is  TIMINT  which  drives  both  the  USA  integrator  for 
the  fluid  and  the  STAGS  integrator  for  the  structure.  TIMINT  places  the  structural  dis¬ 
placements  and  velocities  and  the  total  fluid  pressures  into  mass  storage  at  each  time 
step.  User  selected  response  data  are  then  displayed  in  tabular  and  graphic  form  at  the 
conclusion  of  the  run.  The  processor  POSTPR  has  the  same  display  capability  as  TIMINT, 
but  is  used  for  more  detailed  examination  of  the  results  at  a  later  time. 

In  the  following  these  separate  functional  components  are  discussed  individually. 
Detailed  information  for  utilization  of  the  code  is  contained  in  the  Appendices. 

3.1  STAGS  PREPROCESSING 

Execution  of  a  purely  structural  analysis,  linear  or  nonlinear,  is  carried  out  by 
two  independent  programs  STAGS1  and  STAGS2.  The  STAGSl  program  is  basically  a  prepro¬ 
cessor  which  reads  input  data  defining  a  structural  model  with  loads  and  boundary  con¬ 
ditions  and  prepares  an  output  file,  TAPE2,  for  use  by  STAGS2.  STAGS2  performs  the 
actual  analysis  and  reads  no  input.  The  STAGSl  preprocessor  program  has  been  slightly 
modified  for  use  with  the  USA  Code.  Input  data  for  STAGSl,  as  described  in  the  STAGS 
Users  Manual  [3],  remains  unchanged  except  for  one  logical  record,  B-l.  An  additional 
data  word,  IFLUA,  is  added  to  this  record.  A  value  of  -1  for  this  item  signifies  the 
USA-STAGS  combination  will  be  executed.  The  user  must  also  ensure  that  the  local  co¬ 
ordinate  systems  associated  with  each  wet  surface  of  the  structure  must  be  consistent 
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Organization  of  the  USA-STAGS  Code 


with  normal  displacements  that  are  positive  into  the  fluid.  The  STAGS1  processor  then 
outputs  two  additional  files,  TAPE18  and  TAPE19.  TAPE18  contains  a  list  of  global  co¬ 
ordinates  of  the  structural  nodes.  TAPE19  contains  diagonal  mass  components  for  each 
structural  degree  of  freedom  and  also  a  vector  which  identifies  the  degree  of  freedom 
numbers  with  displacement  components  at  the  nodes.  Both  files  are  sequential  and  may  be 
saved  on  tapes  (by  system  copy  utility  routines)  or  as  permanent  disc  files;  however, 
for  immediate  access  by  USA  processors,  disc  files  are  preferable.  These  data  are  re¬ 
trieved  by  the  fluid  processors  FLUMAS  and  AUGMAT  whose  functions  are  described  below. 
The  STAGS1  preprocessor  also  prepares  a  sequential  "model  definition  file"  on  TAPE21. 
This  file  may  be  used  as  input  with  the  STAGS  plot  program  (STAPL)  to  obtain  plots  of 
the  structural  model. 

3.2  THE  FLUID  MASS  PREPROCESSOR  FLUMAS 

This  code  constructs  the  fluid  mass  matrix  for  a  structure  submerged  in  an  in¬ 
finite,  inviscid,  incompressible  fluid  by  the  boundary  element  technique  [16].  In  addi¬ 
tion,  FLUMAS  can  form  the  mass  matrix  for  a  body  in  the  vicinity  of  a  fluid  free  surface 
through  the  use  of  imaging  techniques.  Out-of-core  processing  is  facilitated  by  use  of 
the  virtual  memory  system  VMSYST  so  that  core  size  is  not  a  limitation  on  rhe  number  of 
fluid  DOF.  The  code  also  generates  fluid  mesh  data  and  a  set  of  transformation  coeffi¬ 
cients  relating  the  structural  and  fluid  DOF.  The  computation  of  these  coefficients  is 
based  upon  the  use  of  centroidal  nodes  for  the  fluid  elements  and  the  assumption  of  a 
bilinear  variation  of  displacement  over  the  surface  of  each  structural  element.  This 
assures  that  the  description  of  the  fluid  pressure  forces  in  the  two  mesh  systems  is 
statically  equivalent  without  inducing  moments  at  the  structural  nodes.  Finally,  the 
code  generates  the  symmetric  matrices  Dj^  and  D^  that  appear  in  the  computational  form 
of  the  DAAj  and  DAA2  equations  which  involve  the  inverse  of  the  fluid  mass  matrix. 

FLUMAS  contains  a  refined  formulation  for  the  fluid  mass  matrix  that  includes  the 
primary  effects  of  element  curvature.  In  addition,  it  has  the  capability  to  treat 
structures  containing  both  general  geometry  and  arbitrary  axis,  multi-branch,  multi¬ 
harmonic  surface-of-revolution  components,  as  described  in  [10],  The  code  can  also 
efficiently  construct  the  fluid  mass  matrix  for  a  body  with  one  or  two  planes  of  sym¬ 
metry  by  using  a  mesh  which  covers  1/2  or  1/4  of  the  surface.  Symmetric  or  anti¬ 
symmetric  fluid  motions  can  then  be  imposed  on  the  portions  of  the  surface  not  covered 
by  the  mesh.  Two-dimensional  "plane-strain"  behavior  of  long  cylinders  can  also  be 
simulated.  The  code  contains  an  automatic  mesh  generator  for  cylindrical  surfaces  and 
an  improved  error  exit  control  that  facilitates  fluid  mesh  debugging.  Finally,  a  useful 
diagnostic  tool  in  the  code  is  the  capability  to  solve  the  fluid-boundary-mode  problem 
Mf  _u_  ■  AAf_u_  [16]. 
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Typical  input  data  for  this  processor  includes 

o  Mesh  geometry 

Fluid  wet-surface  mesh 
Structure  wet-surface  mesh 

o  Element  definitions 

General  curved  surface 
Surface  of  revolution 

o  Material  property 
Mass  density 
Speed  of  sound 
DAA2  parameter 

o  Constraints 

Location  of  free  surface 
Half  model 
Quarter  model 
Long  cylinder 

Node  reassignment  in  fluid-structure  transformation 
A  detailed  description  of  the  required  input  data  is  given  in  Appendix  C. 

3.3  THE  AUGMENTED  MATRIX  PREPROCESSOR  AUGMAT 

This  processor  uses  data  from  both  STAGS1  and  FLUMAS  to  construct  the  specific 
matrices  required  for  solution  of  the  augmented  Eqs.  (2.9)  or  (2.34)/(2.35).  The  output 
of  this  code  includes  not  only  the  required  matrices  in  skyline  form,  but  also  a  conden¬ 
sation  of  the  output  from  both  STAGS1  and  FLUMAS.  This  has  been  done  so  that  only  one 
permanent  file  need  by  referenced  as  input  to  the  time  integrator;  this  results  in  im¬ 
proved  data  handling  and  core  usage.  In  contrast  to  earlier  versions  of  USA-STAGS, 
AUGMAT  does  not  form  the  fluid  matrices  and  but  rather  puts  them  in  the  skyline 
format  required  by  SKYPUL  [21,22].  and  are  now  formed  only  in  FLUMAS.  Input  to 

this  code  involves  the  following  information 


o  Mass  matrices 
Fluid 


Structure 

o  Structural  DOF  correspondence  table 

External  and  internal  node  descriptions 

Factorization  order 

DOF  reduction  due  to  constraints 

o  Fluid  mesh  geometry 

Global  coordinates  of  fluid  nodes 
Direction  cosines  for  nodal  surface  normals 
Areas  of  fluid  element 

o  Fluid/structure  DOF  transformation  coefficients 

o  Fluid  material  properties 
DAA2  parameter 

o  Constraints 

Half  model 
Quarter  model 

Although  this  constitutes  a  substantial  amount  of  information,  almost  all  of  it  is  re¬ 
trieved  from  permanent  data  files.  A  detailed  discussion  of  the  required  input  data  is 
contained  in  Appendix  D. 

3.4  THE  TIME  INTEGRATION  PROCESSOR  TIMINT 

This  main  processor  constitutes  an  implementation  of  the  unconditionally  stable 
staggered  solution  technique  developed  in  [4]  for  DAAj.  The  primary  output  is  a  set  of 
permanent  data  files  that  contain  nodal  values  for  structural  displacement,  structural 
velocity  and  wet-surface  pressure  at  every  time  step.  In  addition,  parallel  files  are 
created  that  retain  restart  information  at  time  intervals  dictated  by  the  user.  The 
code  has  a  variable  time  stdp  capability  and  can  treat  a  spherical  incident  wave  of 
arbitrary  pressure  profile  and  source  location.  Exponentially  decaying  waves  can  also 
be  treated  by  providing  magnitude  and  decay  information.  In  addition,  incident  wave 
pressure  and  particle  velocity  are  tabulated  and  displayed  with  a  "printer-plot" 
package.  If  the  body  is  in  the  vicinity  of  a  free  surface,  unloading  due  to  reflection 
of  the  incident  wave  from  the  surface  is  included  and  the  effects  of  bulk  cavitation  on 
the  free  field  pressure  history  are  approximately  treated.  Finally,  selected  response 
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histories  can  be  listed  and  then  displayed  for  immediate  examination  using  the  “printer- 
plot"  graphics  package  embedded  both  in  T1MINT  and  in  POSTPR  (see  Sec.  3.5). 

The  computational  strategy  for  the  staggered  solution  procedure  is  embodied  in  the 
following  eight  steps,  assuming  the  solution  is  known  at  time  t: 

(1)  Estimate  the  unknown  structural  restoring  force  vector  at  t  +  At  from  the 
extrapolation  of  current  and  past  values 

(2)  Transform  this  extrapolation  into  fluid  node  values  and  form  the  right-hand 
side  of  the  fluid  equation,  which  also  involves  the  known  incident  pressure 
at  t  +  At 

(3)  Solve  the  fluid  equation  and  obtain  a  preliminary  estimate  of  the  total  pres¬ 
sure  vector  at  t  +  At 

(4)  Transform  fluid  pressures  into  structural  nodal  forces 

(5)  Solve  the  structural  equation  for  the  displacement  and  velocity  vectors  at 
t  +  At 

(6)  Transform  the  computed  structural  restoring  forces  at  t  +  At  into  fluid  node 
values  and  reform,  the  right  hand  side  of  the  fluid  equation 

(7)  Re-solve  the  fluid  equation  and  obtain  refined  values  for  the  total  pressures 
at  t  +  At 

(8)  Save  system  responses 

Steps  1,  3,  and  5  constitute  the  basic  staggered  solution  technique,  while  Steps  2  and  4 
are  required  because  of  the  difference  between  the  fluid  and  structural  surface 
meshes.  The  iteration  on  the  fluid  solution  reflected  in  Steps  6  and  7  has  been  added 
to  enhance  accuracy.  Inasmuch  as  the  computation  time  is  overwhelmed  by  the  structural 
solution,  this  requires  only  a  small  increase  in  total  run  time.  The  use  of  a  three- 
point  extrapolation  method  in  Step  1  also  improves  accuracy,  as  discussed  in  [4],  The 
3-step  implicit  Park  [23]  integration  method  is  used  for  the  fluid  solution,  while  the 
Integration  algorithms  embedded  in  STAGS  have  been  described  earlier. 

Typical  input  to  this  processor  includes 

o  Incident  wave  characteristics 
Location  of  source 
Location  of  standoff 
Pressure  profile 
Linear  interpolation 
Cubic  spline,  fit 
Exponentially  decaying  wave 


3-6 


o  Time  step  information 

Start  and  finish  times 
Time  increment  values 

o  Restart  data 

o  Display  directives 
Displacements 
Velocities 
Pressures 

Detailed  user  information  concerening  TIMINT  is  given  in  Appendix  E. 

3.5  THE  RESPONSE  POSTPROCESSOR  POSTPR 

This  utility  is  responsible  for  the  listing  and  "printer-plot"  as  well  as  "vector- 
plot"  graphic  display  of  selected  system  responses  and  pseudo-velocity  shock  spectra. 
Output  files  containing  the  structural  displacement  field  at  user-specified  instants  in 
time  may  also  be  created  from  the  response  history  files  to  provide  "snapshots"  of  the 
deformed  structure.  Some  of  the  same  capabilities  are  also  embedded  in  the  TIMINT  pro¬ 
cessor  for  immediate  selective  scanning  of  the  output.  POSTPR,  however,  is  used  for 
more  detailed  examination  of  the  results  at  a  later  time.  As  a  complete  display  of  all 
structural  and  fluid  DOF  histories  for  even  a  moderate  size  problem  could  run  into 
thousands  of  pages  of  output,  care  must  be  exercised  in  the  selection  of  data  to  be 
displayed.  Usage  of  this  code  is  discussed  in  Appendix  F. 

3.6  THE  DATA  MANAGER  DMGASP 

DMGASP  is  a  self-contained  utility  module  that  functions  as  a  manager  of  auxiliary 
storage  and  as  the  focal  point  for  all  block  input/output  activities  [24],  Constituting 
the  lowest  level  of  the  NOSTRA  Data  Management  System  [25],  it  carries  out  the  direct 
transfer  of  data  blocks  between  core  and  peripheral  storage.  (The  terminology  "direct 
transfer"  is  used  here  to  denote  unformatted  and  unbuffered  data  transmission.)  The 
basic  auxiliary  storage  management  operations  embodied  in  DMGASP  are 

o  Activate  storage  device 
o  Position  device 
o  Read  data  block  from  device 
o  Write  data  block  on  device 
o  Deactivate  device 
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In  the  USA  Code,  DMGASP  Is  operated  as  a  stand-along  I/O  package  that  receives 
directives  directly  from  the  master  processors.  Assembly  language  versions  of  DMGASP 
currently  exist  for  UNIVAC  1100  EXEC-8,  CDC  SCOPE  3. A  (NOS/BE),  rnd  CDC  NOS  operating 
systems. 

The  UNIVAC  version  of  DMGASP  has  been  operational  since  197A  and  embodies  a  wide 
spectrum  of  extensively  tested  functions,  including  magnetic  tape  handling.  On  the 
other  hand,  the  CDC  SCOPE  version  has  been  operational  since  mid-1977  and  has  been 
tested  only  on  a  subset  of  functions  dealing  with  mass  storage  I/O.  The  CDC  NOS  is  a 
variant  of  the  CDC  SCOPE  version  and  has  been  used  since  late  1977. 

Finally,  a  FORTRAN  77  version  of  DMGASP  was  produced  in  1980  for  the  VAX  11/780. 

3.  THE  VIRTUAL  MEMORY  SIMULATOR  VMSYST 

VMSYST  is  a  virtual  storage  simulator  for  computers  that  are  not  built  around  a 
rtual  memory  system  [26].  All  data  in  the  virtual  system  is  partitioned  into  pages, 
ich  are  blocks  of  consecutive  data  words  of  a  fixed  page  size.  Pages  residing  in  core 
the  page  buffer  are  called  active  pages.  Inactive  pages  are  resident  in  auxiliary 
storage  only.  In  this  utility  the  page  and  page  buffer  sizes  can  be  conveniently  ad¬ 
just  1  to  suit  the  application.  Input  and  output  to  auxiliary  storage  is  handled  by 
DMGASP;  otherwise  VMSYST  is  written  in  transportable  FORTRAN. 

The  primary  advantage  of  a  virtual  memory  system  is  the  efficient  processing  of 
many  small  records  such  as  columns  or  rows  of  large  full  matrices  that  can  be  treated  as 
vectors.  In  essence  VMSYST  keep  track  of  whether  a  desired  block  of  data  is  resident  in 
core  in  the  page  buffer,  or,  has  been  moved  to  an  external  storage  device  by  DMGASP.  If 
it  is  not  currently  resident  in  the  page  buffer,  VMSYST  retrieves  it  and  makes  it  avail¬ 
able  to  the  application  progrm.  This  double  movement  of  data  is  the  major  price  paid 
for  the  benefits  of  the  virtual  system. 

In  USA-STAGS,  VMSYST  is  used  for  the  out-of-core  generation  of  the  fluid  mass 
matrix  described  in  Section  3.2. 

3.8  THE  SKY PUL  PROCESSOR 

SKYPUL  (SKYmatrix  Processing  Utility  Library)  is  a  system  of  computer  codes  de¬ 
signed  to  process  large-order  sparse  symmetric  matrices  stored  in  the  so-called  skyline, 
profile,  envelope  or  variable-band  arrangement  [22].  Although  the  fluid  equation  system 
produce'  by  the  family  of  doubly  asymptotic  approximations  is  full,  SKYPUL  can  process 
this  system  in  a  multi-block  out-of-core  mode  so  that  there  is  no  essential  restriction 
on  the  number  of  fluid  degrees  of  freedom  that  can  be  used. 


The  main  virtues  of  SKYPUL  are  the  generality  of  the  representation  and  the  simpli¬ 
city  of  storage  resource  management.  The  latter  advantage  is  reflected  in  the  clean 
implementation  of  the  basic  skymatrix  processing  operations  such  as  factorization  and 
equation  solving.  The  resulting  codes  are  relatively  easier  to  understand  and  maintain 
than  those  based  on  more  sophisticated  sparse  matrix  storage  schemes. 

A  distinguishing  feature  of  SKYPUL  is  that  all  of  the  constituent  programs  comply 
with  the  organizational,  internal  documentation,  and  array  self-description  standards 
set  forth  in  the  General  Description  Manual  of  the  Engineering  Analysis  System  (EASY) 
[21]. 
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SECTION  IV 


EXAMPLE  PROBLEMS 

This  section  is  eventually  intended  to  contain  a  collection  of  sample  problems  for 
USA-STAGS  that  covers  the  general  range  of  application  to  linear  and  non-linear  analysis 
of  structures  with  simple  geometric  shapes,  i.e.,  cylinders  and  spheres.  At  this  time 
only  one  sample  problem  is  included,  that  of  the.  infinite  circular  cylindrical  shell 
excited  by  a  transverse,  plane,  step  wave.  The  case  considered  here  is  for  linear  elas¬ 
tic  behavior  only  and,  although  this  does  not  fully  demonstrate  the  computational  power 
of  USA-STAGS,  it  does  lend  itself  to  direct  comparison  of  the  same  problem  treated  by 
USA  in  its  stand  alone  configuration  [2]. 

4.1  SUBMERGED  INFINITE  CYLINDRICAL  SHELL 

For  this  oroblem,  a  72-node,  36-element  model  with  a  uniform  circumferential  mesh 
was  constructed.  The  length  of  the  cylindrical  shell  equalled  the  circumferential  di¬ 
mension  of  the  square  plate  elements  used  for  the  model;  hence  the  shell  was  one  element 
long.  Kinematic  constraints  of  zero  axial  displacement  and  no  end  rotation  were  en¬ 
forced  through  appropriate  input,  to  STAGS1.  The  fluid  model  consists  of  36  equally- 
spaced  elements  around  the  circumference;  the  two-dimensional  nature  of  the  infinite 
shell  geometry  was  simulated  by  exercising  an  option  in  the  fluid  preprocessor  FLUMAS 
that  adds  fictitious  elements  in  the  axial  direction. 

The  transverse,  plane  step-wave  is  of  unit  incident  pressure  and  material  proper¬ 
ties  are  used  that  correspond  to  a  steel  shell  immersed  in  water.  The  input  data  are 
normalized  so  that  the  density  and  speed  of  sound  for  the  fluid  both  equal  unity;  hence, 
the  density.  Young's  modulus,  and  Poisson's  ratio  for  the  structural  material  are  taken 
as  7.85,  98.125,  and  0.3,  respectively.  The  radius  and  wall  thickness  of  the  cylinder 
are  1  and  0.01,  respectively. 

Computational  results  are  shown  in  detail  in  Appendices  B  through  F  and  it  is  noted 
that  there  is  little  discernable  difference  between  the  USA-STAGS  and  USA-SPAR  analyses 
(See  [2]). 
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REDUCE (- ) 

USAS.  READ  USAS  INPUT 

A-4 


******* 

CCNTRCL  CARC 


STAGS  1“ USAS  EXECUTION  (RESTART) 
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SAMPLE  PROBLEM  FOR  STAGS1  PREPROCESSING 
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This  appendix  contains  a  sample  input  deck  and  subsequent  output 
for  the  infinite  cylindrical  shell  problem  presented  in  Section  4.  See 
STAGS  user  manual  [3]  for  detailed  instructions  for  deck  setup. 
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APPENDIX  C 

USER  INFORMATION  FOR  THE  FLUID  PREPROCESSOR  FLUMAS 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  in  Section  4. 
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REASONABLE  CIRCUMFERENTIAL  SPACING  OF  THE 
ELEMENTS.  THE  LATTER  CAN  GENERATE  A 
SLIGHTLY  MORE  ACCURATE  FLUID  MASS  MATRIX 
HOWEVER  THE  FORMER  CAN  GIVE  A  SLIGHTLY 
BETTER  STRUCTURAL  RESPONSE  CALCULATION 
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The  following  discussion  is  provided  ss  sn  aid  to  user  understanding  of  the  Maple 
output  that  is  included  have* 

The  first  item  needing  explanation  is  the  block  subdivision  table.  During  construe* 
tion  of  the  mass  matrix  the  code  must  determine  whether  a  particular  fluid  DOF  pertains 
to  a  GEN  (includes  both  GEN  and  CYL  elements)  or  SOR  element.  Zn  the  latter  case.  It 
oust  also  store  the  branch  or  axis  of  the  element,  its  harmonic,  and  also  whether  that 
DOF  corresponds  to  a  cosine  or  sine  function.  DOF  with  similar  characteristics  are 
■naturally  kept  together  in  the  some  block.  When  the  mass  matrix  is  automatically  pro* 
cessed  in  an  out  of  core  mode  GEN  elcMnts  are  also  partitioned  into  blocks  for  compute- 
tional  convenience. 


The  paraaatars  appearing  in  the  block  subdivision  table  are: 

ISUB  •  block  number  j 

ITYP  -  GEN  or  SOR  \ 

IBEG  -  first  row  of  block 
IROW  -  number  of  rows  in  block 
1BRA  -  SOR  branch  or  axis 
I  HAP.  -  harmonic  number 

I FUN  -  COS  or  SIN  »  \ 

\  \ 

Next,  the  terms  appearing  under  /'Fluid  Hash  Geometric  Arrays'*  are  defined  as: 


X.Y.2 


NX.NY.NZ 


-  number  of  comer  points  for  a  particular  fluid  element 

-  global  cartesian  coordinates  of  the  fluid  element  centroldsl 
control  point 

-  components  of  the  outward  unit  normal  veetor  for  the  fluid 
element 

*  number  of  structural  node  points  that  are  coupled  to  a  p ar¬ 
ticular  fluid  element  for  the  purpose  of  force  application 


AO0,A20,All,A02  -  area  and  moments  and  product  of  inertia  of  fluid  element* 
Used  Internally  for  construction  of  the  fluid  mass  matrix 
and  of  the  fluid-structure  transformation  coefficients  for 
general  elements.  For  SOR  elements,  these  values  sre  for 
the  sub-elements. 


BII.CIl 


-  diagonal  terms  of  B  and  C  matrices  used  for  the  construction 
of  fluid  mass  matrix  (see  (16]) 


When  SOR  elements  are  Included  in  the  fluid  mesh  the  following  new  terms  will  appear 


In  the  output: 


NSOR  -  number  of  SOR  element 

NFLU  -  DOF  in  fluid  mass  matrix 

RAD  -  radius  of  fluid  element  control  point  from  axis  of  revolution 

NCIR  -  number  of  integration  points  or  sub-elements  used  in  circumferential 
direction 

Local  Fluid-Structure  Transformation  Coefficients  appear  next.  This  is  a  summary 
that  indicates  which  structural  nodes  couple  with  a  particular  fluid  control  point  and 
the  weighting  factor  for  each.  The  weighting  factors  must  always  sum  to  unity  for  any 
fluid  control  point. 

The  generalized  areas  that  follow  are  simply  A00  for  GEN  elements.  For  SOR  elements 
with  IHAR  **  0  they  are  A00*NCIR;  for  all  other  SOR  elements  they  become  . 5*A00*NCIR. 

The  eigenvalues  and  eigenvectors  that  follow  the  listing  of  the  added  mass  matrix 
correspond  to  the  "Fluid  Boundary  Mode"  problem  [16].  For  the  infinite  cylindrical 
shell  problem  presented  here,  the  exact  eigenvalues  should  behave  as  1/n  with  correspond¬ 
ing  modes  cos  n6  and  sin  n0.  The  first  eigenvalue  listed,  0.11831+04,  is  an  approxima¬ 
tion  to  «*>  for  n»0  and  it  can  be  seen  that  the  subsequent  eigenvalues  are  relatively  well 
behaved . 

If  a  table  labeled  "SUMMARY  OF  1-0  ACTIVITY"  appears  in  the  output,  this  indicates 
that  automatic  out-of-core  processing  has  taken  place.  In  such  a  case  the  "Fluid  Boundary 
Mode"  problem  is  not  solved  and  its  diagnostic  characteristics  are  unavailable  to  the  user. 
If  there  are  any  serious  errors  in  the  fluid  mesh  geometry  that  have  remained  undetected 
through  the  generation  of  the  mass  matrix  these  may  show  up  in  the  construction  of  the 
matrix  [see  Eq.  (2.6)],  i.e.,  the  occurrence  of  factorization  errors  for  the  elements 

in  question. 

The  following  input  and  output  for  the  infinite  circular  cylindrical  shell  problem 
contain  some  minor  differences  due  to  the  fact  that  the  input  is  appropriate  to  the 
standard  CDC  or  UNIVAC  USA-STAGS  version  3  whereas  the  output  is  from  the  VAX  virtual 
memory  machine.  The  basic  reason  for  this  is  that  the  VAX  version  does  not  explicitly 
process  the  fluid  equation  system  in  a  multi-block,  out-of-ccre  mode  in  contrast  to 
the  CDC  and  UNIVAC  versions.  In  addition,  permanent  file  naming  conventions  differ 
slightly;  however  it  is  anticipated  that  these  differences  should  not  prove  to  be  a 
difficulty  for  the  user. 
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APPENDIX  D 

USER  INFORMATION  FOR  THE  AUGMENTED  MATRIX  PREPROCESSOR  AUGMAT 


This  appendix  Includes  a  copy  of  the  users  manual,  and  a  sample  Input  deck  and 
subsequent  output  for  the  Infinite  cylindrical  shell  problem  presented  in  Section  4. 


l/i  Uii/l  Q 

3  H  -I  _J  UJ 

X  Of  *  <  O 

o  «  u.  a  uj  -j 

we  5  3  t  a 

H  OUh-Zd 

2  UJ  I/1UUI4 

+*  3  tt  5  5  X 

<  o  o  »  a  o  uj 

Ct  ot  h  3  CD 

h  l/l  U.  UJ  |/l  4 

in  uj  -j  I  h 

ZO«H_ji±Jin— ' 

ogxojiD* 

u  z  z  *-  z  c 

u.  *  O 

J  JnQmBI^UI 

<  <  UJ  >-  Ui  H  CD 

2  a  m  in  <  h  z 

03Z3h2*-U* 
O  O  UJ  <  UJ 

h  o  »-*  o  a  am 

«  3  >-  Z  _  £  ~ 

jocH^zom 

1/IKQUJ02ZO 

z  m  z  co  o 

<  o  moo« 

aouuii-D 

t-  h  a  2  m 

V  ^  H  1/1  1/1  o 

u.  o  a  <iwu 

wWI-WBOUh 
HU JH  DhQ 

u  JXUI/IUBZ 
DaiQZUl- 
tta>Q  o  a  *  u 
H^mxou.x*-* 


UJ  UJ  O. 

CD  X  I  uj  5 

u  h  m  x  uj  i 

o  ►<  uj  a  x  : 

nig  X  «  <  h  * 

J2hi-a 
uj  <  m  in  uj  • 

a  z  a  u.  y  Z  i 

irt<wQMO>Z>l 
uj  a  u.  j  uj  o  i 
u  in  uj  o  xzxX 
m  z  m  o  uj  a  o  mi 
ttO  JUJlO  *  I 
HM<ttuu.uinin 

<  h  u.  x  2  >  ; 

X  <  a  '  uj  h  o  <  i 

Sui  uj  o  m  -h  3  « 

l/IHUtt^K  Jl 
ouiMi/)KUjw<<: 

SoSS~Z*-3mi 

cx  *-«  uj  a  m  h  u  h  i 
o  3  X  x  a 

U.  -J  H-  -•  U  Q  D  >  i 
liOttUJlZHH 

o  »-  tt  H  <  3  »-< 


a  m 

UJ  UJ  < 

<  -j  > 

H  w 

H  I  _l  X 

I/U/IU.  H 

m 

Z  H  -J 

m  *-i  m  h 

z  o 

uj  3  uj  X  <  uj  o 

o 

X  2  u.  X 

a  z  x  x  uj 

V  11 

<3  H 

UJ  S  HZ 

3  Z 

h  o  -j  ^  2 

m 

<  52  ^  > 

x  <  z  o 

m  o 

Z  =J 

J  W  ot  H  O 

u.  *■*  o  z  <  d  m ; 

hkuozxhxi/i: 

<  H  Q  M  O  N  u 
uuZuODOftHi 
3o.Hin<uaH 
oraof«*-‘*Hsu.i 
H«aooxmxo: 


o 

z 

H 

O  UJ 

►-«  m 

m 

H  «J 

<  < 
X  U 

< 

QC 

H 

O  uj 

< 

u  in 

o 

in  »-• 

z  X 

> 

<  a 

a 

a  uj 

H 

H  X 

UJ 

H 

X  uj 

uj  O 

o  m 

a 

UJ  -J 

3  * 

O  « 

H  O 

U. 

U  UJ 

X 

3  H 

m  uj 

tt  in 

uj  m 

H  *-« 

X  H- 

m  -i 

i 

O  a 

O  u; 

HU  UJ 

*-  CD 

3  X 

3 

-J  H 

-1  O 

u  O 

U.  H 

u.  • 

u.  m 

«  Q 

*■*  « 

UJ 

UJ  H 

UJ  < 

3  m 

3  H 

a  H-. 

a  < 

H  -J 

H  Q 

i  <  mo 

<  H 

:  Oc  h  X 
i  uj  o  m 
:  >  z  h  o 


<  a 

U.OI- 

h  < 

m  «  X  z 

<  oc  k  o 

UJ  •— 
Z  a  m  h 
uj  O  <  < 
*  a 

>  <  m  x  o 

HiAaO 
UJO  — 
i  u  2  u>  u 
:  cd  u.  2 

>  u  UJ  D 

i  m  w  X  o 
:  >  h  < 

,  <  1/1  l/l  UJ 

:  *  2 

i  J  J  u  O 

■  <  <  X  -I 

CO  H  < 

•  h  o  _ 
i/i  j  2  g 
:  2  o  •-  z 

I  X  < 


CD  X  UJ 

i  <  uj  uj  X 
I  •-  u  h 

I  Qt  2  H 

:  <  ~  ©  Z 
>  m  2  « 


i  a  -j  o 

*  UJ  Q.  ■ 

2  K  H  ^ 

O  uj  O 

i  **  X  2  U. 

h  <  »-«  O 

<  a  a 

!  X  «  CL  UJ 

ina  3 
iX  -J  -J 

o  m  <  < 

|Q[h|-> 

I  0.  X  2 

I  Q.  U  UJ  < 

>  <  X 

a  <  h 
i  o  o  o  < 
«  »-«  u  Z  X 
)  H  3  H 
j  O  uj  u 
.HO  O 
a  •-»  >  uj 
J  X  o  z  > 

:  >  x  <  a 

■  in  o  uj 

<  >  m 

1  Uj  CD  CD 

j  >-  m  o 

)  -J  *H  2 

jccuujz 
1  3  UI  >  UJ 

j  o  a  «-»  uj 

c  e»  Cl  o  CD 


a  uj  > 

O  CD  < 

-J  u.  X 
-J  Z  -j 
uj  h  m  -a 

3  i/i  U  uj  U 

m  uj  o  o 

CD  H  O  -» 

_J  hh 

<  u.  a  -J 

o  CD  ■  UJ  < 

H-  -J  H  2 

ao  jujo 
UJ  H  UJ  X  O 
I  I  <  < 
a  m  m  a  •- 

m  x  «  o 

uj  j  a 

4  UJ  4  UJ 

m  o  a  X 
a  ~  *  h 
ointt  j 
u  •  Q  <  X 
O  Z  O  H 
>  »-«  in  *-« 

O  u.  -J  3 

<  o  >  m 

a  gna 

3  uj  i  — 

o  3  UJ  H  X 
u  -j  h  m 
<.  ►**  H  2 

;»  Z  <  O 
H  ~  X  ~ 
i/>  <  U.  H  H 

uj  Z  < 

CD  UI  H  2  j 

-J  3  UJ 
uj  *■*  uj  o  a 
X  x  X  X  . 


^-JScvSflSSSrjSSnnrjnnnSnnpijjnjvnnnnnjiirjiinijiiririiJiJiJiDiJipiPiJJJJIJW^^^ 


m  £  z 
►*  D  w  Z 

_i 

Ui 

8 

•  Ui 

X  00 

o 

>  uj  o  <— • 

Z  X  uj  uj  ip 
<  a  x  a  w 
x>  < 

a  Z  in  x 

UJ 

tt  a 

o  z  5 

u  ox 

in 

O 

a 

—> 

< 

tt 

3 

►~ 

X  Z  O  a 

UI  K  UJ  UJ  If*)  tt 

T.  N*  X  tt  ►-  a  * 

r—  a  a  *-«  z  a 

X 

a 

a 

a  >  z 

Z  J  I 

:  3 

X 

O  J 

in  m 

O  a 

o  o 

O 

a  3  a  a 

in 

*n  a  a  o  a 

Ui  U.  I- 

a 

UJ  *J 

*  «/) 

3  o 

Q  •  a  tt 

O  ui 

3 

a  X  a  O  X  K  O 

> 

ft  I  4  U  I 

S 

ui  a 

_i 

UJ  *-t 

x  Q  Q 

u  < 

uj  in  o  o 

ui  Z 

tt 

2  uj  a  a  «(  a 

m 

*-«  a  Q  a  a 

ui  a  q  a  ui 

< 

DC  X 

ui  ui  uj  2  uj 

o  m  in  y- 

in  uj  tt  a  in 

tt  -i  o 

K 

*-  *-  O  a  -j  X 

a  2  in 

-J  X  Vi. 

l  X 

a 

U 

S  ^  5  -  S 

«-•  >  in 

3  o  n  o 

<  *- 

in 

a  i/)  *-•  a  f-  m 

a 

0  3  O 

ui  i-  a 

:  uj  u 

3 

m 

Ui  < 

wtWU 

UJ  tt  K 

a  >  3  i-  < 

a 

a  a  o  a  a 

*“ 

•  co  a 

♦- 

U.  H  ft 

«-i  UJ 

uj  x  a  x 

uj  a  in  o 

3  3  « 

o 

o  in  a  p  o  in  x 

<  Z 

X  OX 

o  z  if 

»  3 

O 

Ola 

tt  r>  in  in 

ttUJh 

CO  a  u  o 

— J  H-  V— 

Ui 

a  2  *-*  ♦-  3 

z  o 

u  a  ►—  a 

*->*-<  u. 

1  K  U 

3 

a  m 

X  uj  > 

o  ~  <  H 

<  o  o 

^  U  tt 

z 

O  a  3  Z  a 

»— 

<  a  < 

=3  o  m 

a 

.  z 

O  a  q  •  £ 

uj  x  a 

OX  tt  z 

>  3  * 

»-  k  a  in  a  a  a 

O  a 

nuis 

-i  o 

a 

ui  in 

Q  UJ  *  ^  J  l/l 

c  a  • 

a  UJ  O  O 

a  a 

-J 

<  ^  «  a  X 

tt  3 

a  tt  <  a  a 

U.  Ui  If 

IS  B  UJ 

m 

uj  O  w 

UJ  ft  <  JU 

in  o 

o  -j  a 

X  u 

> 

O  2  tt  a  Z 

O  a 

*•*  o  a  #i  o 

in  *~ 

1  O  U 

a  «  O 

UJ  UJ  4(  < 

JWHZ 

*  -»  a  x 

^  in  Z 

a 

a  •-»  a  <  Z  a  •-« 

O  O 

X  2  Z  9  Z 

DC  U  < 

z 

O  a 

a  in  •  —  a 

<  o  in  h 

U  3  ~  tt  tt 

tt  Ui 

m 

O  O  Z  >-  O  a 

u  in 

in  O 

ami 

:  o  ui  _i 

»  in  CD  CL 

** 

ui  in  uj 
^  _  3 

it  *-t  (m  h  •  in 
^  x  in 

o  -i a  cm  a 

-1  u  j  o 

V-  O  UJ 
<  UJ  CO 

o; 

a  tt  a  <  > 

a  o  O  Q  *—  o  a 

a  a  ■■ 

ai/il  ui/il 
a  a  O  a  O 

<  a 

:  m 

in 

X  — J 

u  x  •  a  z  x 

<  u. 

ft)  D  uj  O  uj 

X  z 

o 

2  0  3  *-•  a 

O  «  a 

4  K  O  UJ  h  Q 

a  x 

UJ  Ui  Ui 

a 

— 1  a  < 

0  O  ^  O  a 

a  >  z  a 

O  X  X  tt  N 

*-<  in 

u 

o  tt  m  a  3  < 

tt  tt 

Z  <  a  tt  <  a 

u  rj  U  3  X 

Z 

*-H  J  «-4  *-« 

<  2  <  uj 

tt  uj  u  O 

tt  -J  < 

in  <  •-»  o  a  o  a 

a  3  * 

~  Z  a  3  Z  a 

u.  a  < 

:  o  -i  < 

QC  X 

uj  x  in  <  a  a 

-J  *  a 

ui  z  ui  ax 

O  >  X 

UJ 

KaaX>a*-uU 

Cu  a 

O  h  ft  a  «  a 

a  < 

a  ^  a 

O 

3  uj 

uj  a  x  z  < 

in  uj  z 

a  r-«  X  a  *-« 

N 

a  <  3  ►-  a  v-  < 

>  o  m 

tt  o  a  O  O  a 

x  ui  q  a  > 

a 

a  O  X 

a  ooi-o 

Z  a  X  N- 

-i  O  uj  n  z 

m  in  x 

»-  X 

in  a  u  tt  >  <  a 

a  3  a 

o  tt  3  tt 

*-•  rsi 

U  Z  a 

15  UJ  O  «  o  Z 

<  <  a 

in  >  m  x  10  ^ 

ic 

in  *-* 

3  ^  •  Z  X  tt 

tt  3 

O  O  a  tt  o  a 

a  < 

:  i/>  O  < 

Ui 

3  — 

UiQ  Ulh  tt  h 

tt  z 

Q  X  a  X 

U  uj  tt 

a 

<tauuju<03 

a  a  a 

o  o  O  a  o  o 

H-  U 

in  a  u 

o 

a  o  in 

O  O  UJ  <  (J 

a  x  Z  Ui 

ft  1/)  UJ  M  Uj 

O  X  ►- 

¥  »- 

a  o  tt  o  O  a  in 

X  in  < 

o  mo 

<oc 

1  <  UJ 

o 

►“  <*.  Ui 

2  ft  Jft  < 

Ui  M  X 

o  in  uj  X  <n 

-J  ►“  < 

(J  < 

<  3  <  <  m  3 

a  > 

a  a  >  a 

X  2 

;  *  *  g 

z 

in  u  X 

a  u  in  o  u 

a  o  ^  a 

a  o  d  a  ^ 

CO  X 

O  X 

Ottttaaa<a 

a 

4  J  UJ  O  J  iJ 

UJ  «« 

2  O 

tt  •-< 

in  j  Zutt 

lOZO 

UI  UJ 

X 

-J 

a  tt  u  O  < 

m  x  a 

QQ  <  tt  *-  <  tt 
0  0^300 

Ui  3 

OQh 

u 

U  UI  V— 

UJ  <  •  <  UJ 

UI  N  UJ 

u  Z  tt  uj 

u  u  in 

co  in 

uin>a<3ou 

uj  a  a 

a  J  3 

!  h  i  k 

o 

O  a 

O  a  in  a  o  a 

e>  u  r> 

O  O  3  <  O  N 

O  *->  in 

in 

O  Z  a  a  a  a  *m 

co 

a  o  a  -i  o  a 

3  <  3 

i  3  < 

Z  ^ 

a^Hhiuz 

ft  I  Z^' 

N  J  Q  *-U 

X  uj 

X  uj 

a  <  Z  tt  o  a  tt 

<  Z  < 

O  a  o  a  a  o 

a  >  a 

:  -1  <  a 

a 

a  hui 

<  a  in  >  «-• 

<  Ouj 

tt  h  4  Q  O  </> 

tt  a  z 

2  2 

tt  O  n32<q 

Z  •-*  a 

« 

u  3 

Ui 

Ui  w 

J  O  >  ui  ft 

J  O  ft  UJ 

UJHtt  I  4 

uj  a 

5  a 

a  x  in  m  a  2 

O  CL 

i  i  i 

>  <  m  x  a 
a  in  uj  v-  x 

S 

CO  Z 

X  UJ  x 

3  ^  COuJ  2 
a  a  <  in  uj 

Ui  Ui  Ui 

UJ  UI  IL  ft 

CO  1-  UJ  o  >  UJ 

x  a  z  u  *  •  j 

CO  o  u. 

X  H-  ^ 

a 

X  M 

co  a  o  a  a 

xo.-«inaara 

M  Q  UJ 

in  a  o 

O  (N 

D  u.  <  x  «  o 

UHfhlU 

$ 

is* 

X  a  Z  O  UI  X 
a  in  <  a  a  a 

X  tt  Ui  x 

1-  U  ft  1- 

2  <  ui  (,)  2  n 

2  a  19  <  3  u 

3  Z  K 

Z  «  in 

<  *- 

X  in 

3>X303iU> 

ZajliUDOU 

a  in  o 

O  3  < 

to  < 

52 

UJ  _J 

a  -h 

2 

52 

X  < 

K  Jl- 

to 

O  to  on 

•— « 

Z 

UJ  UJ 

>“Z  to 

g  z  <  z 

*4 

O  >  tt 

tt  H 

UJ  UJ 

^lO  W  tt  h 

> 

X  X 

O  UJ 

ZhBQh 

X 

H  J  O 

ODD 

X  X 

>  o  > 

uj  to  uj  <  O  uj  Q  Z 

< 

> 

H  H 

UJ  *-•  to  Q 

*M  <  O  < 

H 

Z  OH 

u.  -1  to 

H  H 

a  a  uj  oc 

cc  a  u  »-  <  a  u.  uj 

X 

< 

Z  H  HQ 

Q  tt  -j  tt 

to  o  UJ  o 

D 

to 

<  *-*  X  < 

UJ  •—  1/1  ♦-  <  z 

X 

g  z 

<  O  Z  2 

2  h  g  h 

H 

H  <  H  UJ 

h  -Oh 

UJ 

to 

g  z 

QI^Q 

i  >  o  O  to  -J  a  uj 

to 

2  H 

tt  _J  H  H, 

UJ  to  to 

< 

H  z  O  z  >  Z  Ui  X  z 

o 

UJ 

Z  H 

2  h  2 

«  Z  X  X  <  to  O  -j 

H 

to 

*■4 

UJ  to  «  -J 

azNz 

O  h  uj  X  uj 

UI  Z  H4 

o 

Q 

H 

3  0  3 

1  2hUKI/)UUJI/)UJ 

2 

H 

Z  to 

CJ  Z  UJ  a  < 

Uj  Q  O 

UJ 

2<tti-H  JZ<  J< 

z 

Q 

Z  to 

OWhC 

i  -J  *-*  <  O  H  to 

UJ 

z 

g  h 

<  D  H  oc 

ouau 

o 

tt  ~  ~  a  « 

O  — J  uj  tt 

2 

g  uj 

DQ  X  CQ 

i  n  >  uj  1  J  <  ui  o 

X 

UJ 

-  z 

O  OC  a  to  D 

o 

Q 

tO  H  3  X  O  X 

a  Q.  o  h 

-J 

X  fi 

H  -1 

u  to  a  2  u  •- 

UJ 

X 

to  UJ 

H  H  Z  H 

n  -j  uj 

z 

H  to  o  tt  X  O  to 

< 

— J 

to  o 

H  <  H 

tt  uj  Ui  uj  O  h  q  3 

—1 

UJ 

to  X 

X  O  (J 

«  *  z 

Z  uj  o  to  o 

O  >  X  z 

CT 

z 

<  z 

to  2 

U  I  X  u 

i  oaiovu-ott  j 

UJ 

-J 

<  UJ 

O  O  O  U  3 

DcZ>0  JZOftUJHZ 

U  tt  o 

D  O 

tt  o 

< 

<  h  a  < 

<  h  <  3  c  a  u. 

UJ 

— J 

UJ  *-  O  oc 

o  o 

<  O  U  H 

H  tt  U 

H 

u 

3  u 

— J 

w  r  2  H 

*  -J  tO  -JO 

o 

2  uj  a 

OC  UJUIK 

•H  -  > 

tt 

►H  uj  <  UJ  o  <  UJ  UJ 

u 

•— 

H  H4 

z  <  a 

2*22 

COZ3U40V)Z 

o 

—  o 

H  O  UJ  l/l  l/l 

*  H  X  -J 

D  H  •  tt  «  o  tt 

IQJ 

D 

u 

^  a  o 

O  _  c 

i  U  ~  to  CO  U  2  o 

3  to 

■-*  to 

O  H 

5  uj  oc  uj 

M  <  Z 

H 

<  UJ  <  O  O  < 

to  I  2  < 

tt 

tt 

3  tt 

^  3  H 

Ul^lU 

»  *  oc  tt  o  x  *-« 

-I  o 

3  o 

O  H4  to 

O  •-«  tt  X  O 

JUJ  O 

u 

H>  ozx<>-«  — 

H 

o 

tt  o 

O  H  to 

Ui  O 

•-HUOW  J  JK*- 

U  u 

-J  u 

UJ  3*3 

LU  _J  H  H 

•  to  X 

3 

H  h  I/U/)  U  T  IO  J  h 

to 

H 

tt  u  z 

Ui  H  O  U 

l  U  O  U.  >-  g  <  3 

►H 

tt  f 

M  "J 

ocau. 

-  Z  H  - 

a 

I  H  UJ  l/IJ  D 

>-  2 

UJ 

to  UJ 

-<  3 

a  <  ui  d: 

£  3  <  CO  2  -J 

UI  Q 

o 

JU.O 

0.0  •  > 

<  Z  H 

H4<JU<4<OHUJUUJ 

UJ  3 

3 

-J  tt  o 

D  Z  ui  2 

1  Z  h  H  -  »  Z  O  m  O 

a  2 

UJ  2 

a  h 

lO  <  to  J 

itiazov  jh«  a 

Z  I  tt 

tt 

-J 

UJ  mJ 

a  h  h 

h  ^  a  h 

OI^Q-J*-J  > 

o 

oc 

CL  O 

■-  to  Z  Z  UJ 

3  H 

tt  tt  H  o  H 

UJ  H  tt  UJ 

O 

< 

tt  « 

a  ^o 

u  o  tt  u 

1  ^(JhUI^OOKUJ 

X  u 

O  U 

<  H  H 

uj  uj  uj  O  O  U 

«..*  H 

< 

10HHZHI-. 

X  UJ  u. 

X 

> 

5  > 

<  H 

3  a  3 

i  h  «  <  Z  uj  ct 

o 

X  O 

oc 

to  ®  uj  *-»  < 

t  uj  O  U 

h  i/i  <  uj  U  Z 

UJ  O  H  X 

tt) 

X 

O  tt 

a  o  u.  a 

as  JZh  i/i  to 

a 

UJ  to  < 

KHQU. 

>  I  UJ  UJ 

H 

X  uj  x  X  O  D 

Jl-ftD 

z 

tt 

UJ 

UJ 

tt  H  < 

H  O  O  H 

O  <  to  h  Z  u 

O  ui 

oc  ut 

£0  O  K 

<OOHUJtt 

H  H4  tt 

< 

h  a  uj  a 

UJ  <  u 

5 

O  X 

tt  X 

CO  H 

l^u  V 

*  ttUJ  JDmO>0 

3 

o  3 

O  to 

HHIUQhD 

UJ  -j  ^ 

O  -j  H  Ui 

tt  q  tt 

< 

o  < 

_  z  to 

1  LU  1 

UJ  o  u  a  2  *-•  -J 

UJ  -J 

-J 

Oh3 

<  O  Z  -J  to 

x  tt  a  a  uj 

tt  >  uj  uj  to  X 

0<0h 

-J 

tt) 

to 

to 

o  o  z 

O  JWC 

1  X  tt  Q  _I-H  JW 

2  < 

UJ  < 

h  Q 

O  to  DO. 

h  uj  a 

_J 

UJ  -J  H  H 

t-<  -J  u 

< 

z 

tt) 

H  U 

~  <  a  - 

*-•  X  >•  O  3  <  <  U 

o  > 

z  > 

2  X 

H  J  u  a  H 

X  <  UJ 

CO  >  tt  z  O  tt 

D  O  < 

*■* 

o 

UJ 

Z  tt 

«  x 

DUOS 

1  H  >  -1  O  X  u  < 

o 

h  o 

tt  Z  <  «3  UJ 

UJ  H  H 

< 

UJ  4*  UJ  H*  <  H 

-J  u  < 

X 

X 

o  x 

H  o 

J  O  UJ  - 

1  J  *  UJ  U.  DC  M  U. 

tt  UJ 

UJ 

Zu.DC 

O  *"  tt  >  X 

>  UJ  UJ  < 

* 

> JI3 J< 

tt  *“•  • 

< 

tt 

H 

H 

Z  U.  tt 

u.  -J  O  u 

i/i  z  i/i  i>  o  i/i  o  u  a 

O  X 

tt  X 

UJ  O  tt 

<  3  Ct  uj 

«  X  CO  z 

o 

ODUJ3I 

O  Ui  uj 

o 

tt 

UJ  O  tt 

ho  uj  <  uj  u  <  3 

« 

o  < 

X 

tt  oc  H  H  CO  UJ 

x  *  ~ 

— J 

X  U  D  u.  U  H 

Ui  2  -J  tt  Z 

52 

o  o 

X 

i 

X  QC  I  UJ  h  l/l  I  |/l 

H  l/) 

to 

UJ  UJ  UJ 

uj  h  u  uj  X 

O  u  — J 

►H  <  H4 

X  UJ  tt  H 

< 

H 

z 

z 

UJ  tt  tt 

KOuiHOttZO 

to 

H 

a  3  o 

eg  to  3  x  o  h 

U  Z  H  tt 

< 

UJ  H  <HXHO.XO 

to 

H4 

H  H 

a  3  g 

r> 

uihja<uj 

tt  UJ 

tO  UI 

U  JZ 

X  Z  tt  X  UJ 

JO  o 

saw  tt  u 

tt  <  o 

H 

tt 

> 

to  > 

y  -1  Z 

i-i/)a-j<<tt3)-j 

H4  X 

<  X 

z  <  < 

2  o  h  >-  uj  z 

«  C.  to  o 

to 

H<IU.OtLUJ)(UJ 

o 

UH 

< 

<  < 

Z<< 

«3t03>UH<< 

tt  H 

-J  H 

«-•  >  oc 

Z  u  to  IO  z  o 

*  D  o 

*-« 

HO.HH-ttttOOQ.UJtt 

z 

tt 

X 

-J  X 

►H  >  a 

m*Tint0H«oo)O*'Ojr>'<rtntpH«DO>o,~c*o*irirH«c>o**cjo<intpH«oo>o,'fccjro*in<i>HflooiO*-c'tn^in®Hcoo)Q 

nr>n«*>r>C}m*<n*'*<n<nntninintnininininintnipgt0ti>tPtpt0ti>gteHr-HHHHHh'HH««ottcoco«®«®«oG) 

NWWWWNMWWNNCIPIMWMNNWNNWWNWNNPlNNNWNNNNNPIPIWMWNNNNWPINNNriWMNWNril 


D-6 


k-  ! 

;,i 


s. 1 

fj 


S  '1 
:  ’  i 

'] 

*  1 


UJ 

o 

z 

> 

> 

♦- 

UJ 

a.  o 

uj  o 

♦- 

* 

« 

< 

< 

UJ 

< 

— J  u 

z 

>-  a 

UJ 

< 

X  ~ 

< 

in 

CD 

z 

-J 

I 

_J 

3  •-« 

fr- 

fr-  UJ 

00 

fr-  •  U. 

UJ 

UJ  > 

* 

* 

+ 

o 

— • 

UJ 

Cl 

fr- 

CL 

<  h- 

1  N 

UJ 

X  »-* 

Z 

00  UJ  < 

c> 

UJ 

Z  w 

ex 

in 

• 

in 

u.  in 

in 

CL 

o 

in 

K  w  u 

o 

in  3  o 

c 

* 

* 

z 

CL 

4-4 

U.  x 

1—4 

ui  o 

<3 

<  UJ  o 

fr- 

3 

a  a  uj 

in 

o  3  a  a 

_j 

U 

u. 

3  Z 

Q 

o  -< 

a 

Q  Z 

X  X  ♦— 

<  »-  a 

>-  in  uj 

00 

« 

i 

HH  Q 

-J  O  UJ 

ex 

_» 

UJ 

O 

in 

fr- 

O 

in 

*-  <  in 

a 

in  *-  n 

fr- 

>-  ex 

O  fr-  UJ 

< 

in  fr- 

< 

<  •  I 

>  < 

z 

<  K 

a 

z 

in  z 

o 

o  z  O  i 

u. 

in 

* 

* 

in  < 

O  fr- 

a  < 

7. 

X  h- 

DO  *-» 

UJ 

in 

*  UJ 

o 

o 

< 

u. 

ex  o  z  o 

o 

* 

3  O 

O  UJ 

UJ 

uj  x 

a 

tu  t-1 

O 

-J  3 

uj  o  in 

o 

►-H 

-J  UJ 

OH  JO 

z 

♦ 

* 

3 

—  ui  Z 

L> 

fr~ 

o 

o  a  > 

a 

f  f 

O 

a  x 

UJ 

fr“ 

3  X  Z 

O 

U  fr-  -j  z 

fr- 

C5 

u 

UJ 

O  fr*  h 

3 

Z  i n 

z 

3  H  j 

uj  a 

U  X 

3  in  uj 

ex 

1—4 

"I  1-  UJ 

lJ 

UJ  H  H 

UJ 

< 

* 

* 

>-  ex 

«N  u  z 

Q 

*i  m 

O  <  Z 

fr-  □ 

a 

a 

o  in 

o 

*  I 

f- 

ex  o  *  < 

in 

►-< 

UJ 

X  - 

O 

Q  uj 

a 

O  X  O  Z  U. 

CL 

o 

o  o  *-> 

UJ 

z 

U.  5 

Z 

2 

O 

♦ 

* 

O  fr- 

ex  uj 

tt 

a  z 

UJ 

a 

l-H 

z 

ex  ♦-  3 

3 

o 

UJ  O 

IM 

uj  O  uj  uj 

O 

o 

H4  Z 

>  »-  -j 

O. 

u. 

o 

a  in  o 

a  > 

-J 

UJ 

a  ex 

— J 

(J 

O  > 

a 

3  o  o  in 

fH 

u> 

• 

» 

ex  uj 

♦-  in  »-i 

— * >•  i 

Z  a. 

z 

in  z 

a  -i 

-J 

in  uj 

< 

OO-i 

CL 

— J  O  3  U. 

X 

Z  uj  i 

O  1 

-J 

o  ^ 

3 

J  UJ  < 

z 

< 

UJ 

-j  uj  X  in 

> 

-J 

u  z  z 

<  JU 

1—4 

o 

« 

* 

UJ  UJ 

Uitt> 

O  1 

-J 

M  h- 

-J  Z 

uj  O 

z 

CO 

JDh2 

< 

~  o 

UJ 

>  < 

u 

in 

mu 

ex  i 

z 

•— i 

K  IP 

o 

^  u.  O  CD 

in 

*-•  -J  O  O 

X 

z 

uj  Z 

CO 

Z  UJ  X 

UJ 

* 

« 

►- 

fr-  >  Z 

a  i 

< 

* 

< 

a 

*  u  a 

UJ 

z 

i  <  *- 

•— « 

ex 

X  Z  UJ 

x  a  x  ~ 

a 

a: 

o 

*- 

-J  o 

UJ  1 

z 

U  -J 

UJ 

*“•  UJ 

-J  3 

>- 

< 

>  -  fr- 

cx 

o 

fr-  -«  3 

o 

•m  o  i-  ex 

in 

o 

4 

* 

in  >  Z  fr-  •-« 

CX  1 

o 

/— S 

o  < 

M 

-‘HN 

-J  -J 

ex 

o 

/—V  UJ  *-H 

K 

z 

U  3 

a  Z  t- 

u 

ex 

3  a 

-•  Z  fr- 

a  » 

UI 

.j  a 

*-  in 

•-«  < 

UJ 

-  x  Z  O 

< 

O  UJ  < 

z 

H  O  < 

< 

* 

4 

X  3 

UJ  < 

• 

O 

w 

3 

o 

W  o  3  > 

> 

u- 

W  «  O  z 

S 

a 

Z  CD  > 

in 

o 

<  ex  z  z 

o 

< 

O  O  a  a 

Z 

z  • 

-J  fr- 

fr- 

O  1- 

3 

a  -j  o 

UJ 

< 

f-H 

in 

X  uj  < 

z 

UI 

♦ 

* 

<  u 

uj  a  uj 

o 

<  • 

••J 

<  u 

UJ  UJ 

m  o 

ex 

a 

UJ  v-  <  u 

fr- 

D 

uj  O 

< 

O  UJ 

UJ 

♦- 

o 

fr-  o 

►h  3  a 

M 

ex  i 

Z  3 

fr- 

Z  Z  ♦- 

x  cc 

O  fr- 

z  < 

in 

Z 

o  x  ex 

in 

UJ 

fr-  Z  O  X 

X 

Z 

• 

* 

< 

u  u  o 

fr- 

C5  i 

c* 

o  a 

UJ 

Ohuj 

Q t  uj 

u. 

3 

O  X  >  -J 

ft 

3 

a  fr—  UJ 

y 

DC 

m  3  a  »- 

3t 

>— * 

O  z 

*— * 

< 

O  i 

O  fr- 

in 

-j  in 

UJ  N 

O 

fr-  < 

*-4 

UJ  fN 

o 

<  UJ 

QC 

« 

* 

< 

fr-  in  u 

a 

ex  i 

u_ 

<K  in 

in 

li.  > 

fr-  1 

ex 

U.  LL  H  Z 

U. 

N  fr—  1 

o 

O 

-J  •  N  U. 

> 

a 

u  u 

1/1  H  < 

UJ 

a  i 

O 

*-4 

UJ 

Z 

O  *  uj 

z 

o 

u 

o  o  >  a 

O 

<  z 

3 

*-* 

O  o 

3 

« 

•— i 

3  > 

a 

i 

O  O 

CO 

o 

in  in 

-J  o 

o 

~  o 

u. 

UJ 

u.  O 

00 

fr- 

li.  UJ  u. 

z 

UJ 

< 

a  X  3  a  4- 

o 

Z  1 

LU 

UJ 

*— 1 

»M 

UJ  *-* 

<  z 

in 

UJ  >  fr-  Z 

o 

> 

o  >  z 

in 

o  >  o  o 

o 

CD 

* 

• 

UJ  u 

»u  Z 

«  1 

2 

UJ  Z 

fr- 

DW  J 

z 

z 

K 

3  <  U 

< 

< 

u. 

o 

<  z 

z 

X  *-« 

fr-  *•<  3 

o 

<  fr 

— 1 

I  ~ 

a 

t— 1 

JItt 

o  < 

o 

2 

j  juq; 

X 

-J 

UJ  -J  < 

O 

z 

X  -J  UJ  UJ 

UJ 

O  O 

* 

* 

3  X 

u  u 

o 

s  > 

< 

fr-  -J 

UJ 

a 

<  h  111  C 

in  3 

<  a  z  uj 

IL1 

a 

3  CL 

o 

LU  CL  3 

3 

fr- 

in 

< 

2  J  UJ  w  K 

u 

UJ 

a 

> 

> 

I 

z 

>  I  <  UJ 

< 

O 

>  in  z  o 

o 

in 

-j  in  uj 

fr¬ 

< 

O  in  -j  uj 

-J 

< 

♦ 

« 

< 

J  JO 

1 

— J 

>- 

u.  * 

*- 

o 

U  fr— 

hh  m 

UJ 

X 

~  o  Z 

z 

<  •-.  I/O 

ill 

Z  ^  <  X 

< 

in  uj 

o 

x  m 

<  u 

a 

Z  fr 

fr- 

z 

u 

< 

o  m 

o 

V 

<  o  o 

O  3 

a 

< 

<003 

o 

>03 

in 

o 

-•  O  >  fr- 

> 

*-* 

a 

* 

4 

a  o  ui  3 

o 

O  1 

< 

Z 

4 

* 

< 

uj  a 

h-  I 

3 

*-*  fr—  in  uj 

a 

t-H  1 

u 

— i 

C5 

u. 

* 

4 

fr- 

u  in  3  ►- 

UJ 

Z  fr 

*-4 

u 

3 

fr- 

fr- 

a 

3  -J  U 

in 

H  1 

a 

< 

in 

UJ 

4 

4 

UJ 

z  z  a  < 

3 

U  1 

UJ 

fr- 

»- 

a 

3 

u 

x  a 

Ui  1 

z 

a 

a 

u 

* 

4 

X 

3  <  UJ  < 

UI 

Q  1 

3 

a 

a 

Z 

a 

UJ 

-J  t—  X  X 

m 

1 

Z 

t-H 

•-« 

1-4 

4 

4 

o  <  <  u 

Z  i 

< 

z 

< 

U  O  2 

> 

UJ  1 

X 

Z 

z 

♦ 

4 

fr- 

0— » 

< 

_J  1 

a 

o 

a 

M 

< 

*— *  U  UI  CO 

z 

CO  1 

-j 

u 

»- 

a 

4 

♦ 

o 

CO  ^j  — 

O  i 

< 

z 

fr- 

u 

w  a  »m  w 

in 

a  i 

> 

a 

m 

4 

4 

fr- 

UJ  u. 

a 

a  i 

z 

in 

a 

z 

3  fr-  S  Z 

UJ 

i 

z 

4 

* 

a 

X  3  ■  uj 

fr- 

-j  i 

3 

fr- 

-J 

Z 

O  Z  in  uj 

U 

<  i 

-l  X  fr- 

»- 

z 

< 

< 

4 

4 

<  O  fr- 

< 

a  i 

o 

< 

in 

X 

a 

> 

> 

is 

Ull  JI 

a 

UI  1 

O  Z  fr- 

o 

CN 

a 

j 

fr- 

fr“ 

a 

a 

a 

4 

4 

— J 

a  uj  in 

< 

Z  > 

a 

fr- 

< 

K 

■j 

a 

a 

< 

> 

> 

JZ  Jnh 

X 

UJ  l 

IN 

a 

a 

< 

in 

i 

a 

a 

s 

z 

4 

* 

< 

H  <  Li.  UJ 

o 

O  » 

in  u 

a 

o  z 

1 

«-csi»><inu)r^ooo>Q*-Nri^in 

0»0)0)OMnO)0)0»0)000000 

MMNMMisrJNNnnnnnn 


•^wn^iniDr'tooo-'-ojrsrfin^ 
nnnnnnnnnnn 


r'-coo>O*-cNonini£Ncoo>O~«n<*iniDr-toao,*,*wonini0r-co 

'•’--MNWPiNNMNNNnnnnnnnnnn^^tn'Mt'ttn 

nnnfonn^wnwnnnnwnnwnrtnnnnnnwnDnnpi 


The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here . 

After  a  summary  of  the  fluid  mesh  geometry  arrays  (see  Appendix  C)  the  first  item 
needing  explanation  is  that  entitled  "Fluid  Element  Wetted  Freedom  Indicator".  This  is 
simply  a  listing  of  the  input  variable  NDICOS  (see  user  manual)  for  each  fluid  element. 

The  section  "Structural  Grid  Point  Numbers  Associated  With  Internal  Sequence  Num¬ 
bers"  contains  a  correspondence  table  that  relates  the  internal  sequence  numbers 
assigned  by  the  fluid  mass  processor  with  the  external  structural  node  number  assigned 
by  the  user. 

The  next  item  entitled  "Grid  Point  and  Freedom  Number  for  Each  Row  of  Stiffness 
Matrix"  identifies  an  integer  vector  that  is  constructed  during  the  STAGS1  preprocessing. 
For  each  structural  equation  the  entry  in  the  vector  consists  of  ten  times  the  struc¬ 
tural  node  number  plus  the  local  degree  of  freedom  number. 

The  last  item  requiring  explanation  is  the  "Freedom/ Equation  Correspondence  Table". 
This  is  an  integer  matrix  of  6  rows  and  as  many  columns  as  there  are  structural  node 
points.  Any  particular  row  corresponds  to  a  local  degree  of  freedom  number  while  a 
column  corresponds  to  the  internal  sequence  number  for  a  particular  external  node  number. 
The  matrix  entry  for  any  particular  set  of  row  and  column  is  the  structural  equation 
number  for  that  pair. 

Depending  upon  user  input  the  various  fluid  matrices  may  then  be  displayed.  The 
matrix  called  TMIT  corresponds  to  £  [see  Eq .  (2.6)],  while  DFDS  denotes  the  sum  of  JD^ 

and  D  [see  Eq.  (2.6)].  In  DAA„  runs  D,.,  is  labeled  DAAl  while  D,„  is  labeled  DAA2. 

~fl  1  2  ~fl  ~f  2 

The  following  input  and  output  for  the  infinite  circular  cylindrical  shell  problem 
contain  some  minor  differences  due  to  the  fact  that  the  input  is  appropriate  to  the 
standard  CDC  or  UNIVAC  USA- STAGS  version  3  whereas  the  output  is  from  the  VAX  virtual 
memory  machine.  The  basic  reason  for  this  is  that  the  VAX  version  does  not  explicitly 
process  the  fluid  equation  system  in  a  nulti-block,  out-of-core  mode  in  contrast  to  the 
CDC  and  UNIVAC  versions.  In  addition,  permanent  file  naming  conventions  differ  slightly; 
however  It  is  anticipated  that  these  differences  should  not  prove  to  be  a  difficulty 
for  the  user. 
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APPENDIX  E 

USER  INFORMATION  FOR  THE  TIKE  INTEGRATION  PROCESSOR  TIMINT 


This  appendix  includes  a  copy  of  Che  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  in  Section  4. 
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The  following  discussion  Is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

First,  the  amount  of  storage  required  for  the  run  given  in  the  output  refers  solely 
to  the  blank  common  chat  is  set  in  the  main  program,  USAS.  An  error  exit  is  taken  if 
insufficient  storage  is  available  and  the  user  must  see  that  more  is  provided  either  by 
a  recompilation  on  UNIVAC  1100-OS  or  by  a  field  length  request  on  CDC. 

Sector  address  information  for  the  response  and  restart  files  that  is  listed  at 
various  places  in  the  output  is  extremely  important  for  subsequent  restart  runs. 

The  next  item  needing  discussion  is  the  transient  response  tabular  listings.  The 
desired  responses  are  displayed  in  matrix  form  so  that  each  row  contains  the  entire 
history  of  a  particular  degree  of  freedom  except  for  the  first  row  which  is  time.  Each 
column  therefore  contains  the  instantaneous  values  of  the  complete  set  of  response  var¬ 
iables  desired  at  a  particular  time.  Each  row  is  identified  by  the  structural  or  fluid 
node  and  its  degree  of  freedom.  The  letters  D,  V,  and  P  stand  for  displacement,  velocity 
and  pressure,  respectively. 

Although  printer  plots  of  the  transient  response  results  can  be  displayed  as  part 
of  the  run  such  output  has  been  deferred  to  the  post-processing  phase  in  Appendix  F  for 
this  sample  problem. 

The  following  input  and  output  for  the  infinite  circular  cylindrical  shell  problem 
contain  some  minor  differences  due  to  the  fact  that  the  input  is  appropriate  to  the 
standard  CDC  or  UNIVAC  USA- STAGS  veision  3  whereas  the  output  is  from  the  VAX  virtual 
memory  machine.  The  basic  reason  for  this  is  that  the  VAX  version  does  not  explicitly 
process  the  fluid  equation  system  in  a  multi-block,  out-of-core  mode  in  contrast  to  the 
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APPENDIX  F 

USER  INFORMATION  FOR  THE  POSTPROCESSOR  POSTPR 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinte  cylindrical  shell  problem  presented  in  Section  4. 
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XVPNAM  A  NAMES  OF  PESPONSE  FILES  THAI  MAKE  UP  A 

CONTINUOUS  SET  OF  TRANSIENT  DATA,  ORDERED 
CHRONOLOGICALLY 
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The  following  discussion  is  provided  as  an  aid  Co  user  understanding  of  the  sample 
output  that  is  included  here. 

The  input  deck  shown  on  the  next  page  requests  vector  plots  for  both  the  transient 
response  histories  and  pseudo-velocity  shock  spectra.  This  is  appropriate  if  the  D1SSPLA 
plot  package  is  available  at  the  users  installation.  Otherwise  appropriate  modifications 
must  either  be  made  to  use  a  different  plot  package  or  the  input  deck  should  be  modified. 
In  any  case  the  printer  plot  package  is  resident  in  USA  -STAGS. 

The  format  used  for  listing  the  pseudo-velocity  shock  spectra  is  similar  to  that 
used  for  the  display  of  the  transient  response  histories  shown  in  Appendix  E  except  that 
the  first  row  is  now  frequency  rather  than  time. 

The  following  input  and  output  for  the  infinite  circular  cylindrical  shell  problem 
contain  some  minor  differences  due  to  the  fact  that  the  input  is  appropriate  to  the 
standard  CDC  or  UNIVAC  USA-STAGS  version  3  whereas  the  output  is  from  the  VAX  virtual 
memory  machine.  The  basic  reason  for  this  is  that  the  VAX  version  does  not  explicitly 
process  the  fluid  equation  system  in  a  multi-block,  out-of-core  mode  in  contrast  to 
the  CDC  and  UNIVAC  versions.  In  addition,  permanent  file  naming  conventions  differ 
slightly;  however  it  is  anticipated  that  these  differences  should  not  prove  to  be  a 
difficulty  for  the  user. 
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ATTN:  W.  Baker 
ATTN:  A.  Wenzel 

RiO  Associates 

ATTN:  P.  Haas 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 

SRI  International 

ATTN:  G.  Abrahamson 
ATTN:  A.  Florence 
ATTN:  W.  Wilkinson 

Systems,  Science  A  Software,  Inc 
ATTN:  D.  Grlne 
ATTN:  T.  McKinley 
ATTN:  R.  Sedgewlck 
ATTN:  Library 
ATTN:  K.  Pyatt 
ATTN:  T.  klney 
ATTN:  T.  Cherry 

Teledyne  Brown  Engineering 
ATTN:  J.  Ravenscraft 

Tetra  Tech,  Inc 

ATTN :  L .  Hwang 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Defense  A  Space  Sys  Group 
ATTN:  A.  Feldman 
ATTN:  P.  Bhutta 
ATTN:  B.  Sussholtz 
ATTN:  P.  Jortner 
ATTN:  A.  Narevsky 
ATTN:  W.  Llpner 

ATTN:  Technical  Information  Center 
TRW  Defense  A  Space  Sys  Group 


ATTN: 

E. 

Wong 

ATTN: 

F. 

Pleper 

ATTN: 

P. 

Dal 

ATTN: 

G. 

Hulcher 

Weidllnger  Assoc,  Consulting  Engrg 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 

Weidllnger  Assoc,  Consulting  Engrg 
ATTN:  J.  Isenberg 
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